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ABSTRACT 


A large  number  of  the  Chapman  and  Cowling  collision  integrals 
have  been  calculated  for  gases  obeying  a modified  Buckingham  potential, 
cP(r)*l_^[oCe  "Ir/J-  The  results  are  tabulated  over  a large 

temperature  range,  kT/e  from  0 to  200,  and  for  four  values  of  the 
parameter  or  , 12,  13,  14,  and  15.  The  treatment  was  entirely  classical, 
and  no  corrections  for  quantum  effects  were  made.  The  results  should  be 
applicable  to  most  simple,  non-polar  gases,  and  may  be  used  to  obtain  informa- 
tion about  intermolecular  forces  from  the  observed  temperature  dependence  of 
gaseous  transport  properties. 

The  second  approximation  to  the  thermal  diffusion  ratio  and  the 
third  approximation  to  the  coefficient  of  ordinary  diffusion  have  been  derived 
according  to  the  method  of  Chapman  and  Cowling. 

Evaluation  of  the  potential  parameters  for  specific  substances  and 
applications  of  the  results  will  be  published  later. 


* This  work  was  supported  in  part  by  Contract  N7onr-28511  with  the  Office 
of  Naval  Research. 
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I.  INTRODUCTION 

The  influence  of  intermole cular  forces  on  the  transport  properties 
of  gases  is  well  known,  and  information  concerning  intermolecular  forces  may 
be  obtained  from  the  temperature  variation  of  gaseous  transport  coefficients 
by  means  of  the  theory  developed  by  Enskog^  and  by  Chapman^.  Such 
information  is  useful  in  correlating  other  properties  of  gases,  as  well  as 
properties  of  the  liquid  and  solid  states. 

The  theory  of  Enskog  and  Chapman  depends  on  the  following 
assumptions;  (1)  only  binary  collisions  between  molecules  are  important; 

(2)  the  binary  collisions  are  elastic;  (3)  the  intermolecular  force  field  is 
spherically  symmetric;  (4)  molecular  collisions  are  adequately  described  by 
classical  mechanics.  Thus  the  theory  applies  strictly  only  to  monatomic  gases 
at  moderate  pressures  and  at  temperatures  high  enough  that  quantum  effects 
are  negligible.  The  extent  to  which  the  theoretical  relations  do  not  nply  when 
conditions  (2)  and  (3)  are  not  met  is  not  yet  known^.  However,  the  theory  has 
had  success  in  correlating  transport  phenomena  in  polyatomic  gases,  so  that 
these  conditions  may  not  be  a severe  limitation^. 

^ D.  Enskog,  Phys.  Zeit.  12,  56,  533  (1911);  Inaug.  Diss.  Upsala  (1917). 

^ S.  Chapman,  Phil.  Trans.  Roy.  Soc . (London)  A211,  433  ( 1912);  A216, 

279  (1916);  A217,  115(1917). 

/ 3) 

' ' In  the  case  of  the  thermal  conductivity  it  is  known  that  the  Enskog- Chapman 
theory  fails  badly'  when  energy  may  be  transported  by  means  of  the  molecular 
internal  degrees  of  freedom. 

S.  Chapman  and  T.  G.  Cowling,  The  Mathematical  Theory  of  Non-Uniform 
Gases  (Cambridge  University  Press,  Teddington,  England,  2nd  edition,  1952), 
p.  7. 


(4) 
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In  order  to  obtain  specific  information  about  inter  molecular  forces 
from  transport  phenomena,  it  is  necessary  to  assume  an  analytical  form  for 
the  intermolecular  potential,  calculate  from  this  the  Enskog-Chapman 
"collision  integrals'1,  or  "temperature. -dependent  cross  sections",  as  a 
function  of  temperature  and  the  potential  parameters,  and  evaluate  the 
parameters  for  any  specified  gas  by  comparison  with  collision  integrals 
derived  from  experimental  measurements.  If  the  theory  is  correct  and  the 
chosen  potential  form  suitable,  a single  set  of  potential  parameters  should 
suffice  to  describe  all  the  transport  properties  of  the  gas,  as  well  as  other 
properties,  such  as  the  equation  of  state.  Except  for  very  simple  and  physically 
unrealistic  potential  forms,  the  calculation  of  the  collision  integrals  involves 
extensive  numerical  integrations. 

The  best  intermolecular  potential  used  to  date  for  the  study  of 
transport  phenomena  is  of  the  Leonard- Jones  form,  with  a repulsion  term 
varying  as  the  inverse  twelfth  power  of  the  distance  of  separation  between  centers 
of  two  molecules,  and  an  attraction  term  varying  as  the  inverse  sixth  power  of 
the  separation  distance.  The  collision  integrals  for  this  potential  have  been 
evaluated  independently  by  a number  of  workers^-^8),  and  their  results  have 

(5) 

T.  Kihara  and  M.  Kotani,  Proc.  Phys.-Math.  Soc.  Japan  j!5,  602  (1943). 

^ J.  de  Boer  and  J.  van  Kranendonk,  Physica  14,  442  (1?48). 

W Hir8chfeider,  Bird,  and  Spotz,  J.  Chem.  Phys.  16,  968  (1948); 

1343  (1949). 

J.  S.  Rowlinson,  J.  Chem.  Fhys.  17,  101  (1949). 
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(9)-(14) 

been  extensively  compared  with  experiment  . The  results  obtained 

have  not  been  entirely  satisfactory,  and  the  question  arises  whether  this 
is  primarily  the  fault  of  the  theory  or  of  an  inadequate  potential  form.  That 
the  Lennard- Jones  (12-6)  potential  is  not  entirely  adequate  is  shown  by  lack 
of  complete  agreement  between  theory  and  experiment  even  for  the  rare  gases. 

It  is  the  purpose  of  the  present  article  to  describe  calculations  of 
the  Enskog -Chapman  collision  integrals  for  a more  realistic  potential  form 
than  has  been  used  previously.  It  is  hoped  that  some  light  mav  be  thrown  on 
the  reasons  for  the  present  lack  of  complete  agreement  between  theory  and 
experiment.  The  inverse  sixth  power  term  in  the  Lennard- Jones  (12-6) 
potential  represents  the  leading  term  in  the  theoretical  form  for  the  dispersion 
energy,  but  the  inverse  twelfth  power  is  only  an  empirical  approximation  to 
the  repulsion  energy,  and  is  not  expected  to  be  accurate  over  a large  range  of 
separation  distance.  There  is  a considerable  amount  of  theoretical  and 
experimental  evidence  which  indicates  that  the  repulsion  energy  is  more 


(9)a  jjirschfelder,  Bird,  and  Spotz,  Chem.  Rev.  44,  205  (1949);  Trans.  Am. 
Soc.  Mech.  Engrs.  J7_l,  921  (1949). 

k Hirschfelder,  Curtiss,  Bird,  and  Spotz,  The  Molecular  Theory  of  Gapes 
and  Liquids  (John  Wiley  and  Sous,  New  York,  1954),  Chapter  8. 

K.  E.  Grew,  Proc.  Phys.  Soc.  (London)  62,  655  (1949);  K.  E.  Grew  and 
T.  L.  Ibbs,  Thermal  Diffusion  in  Gases  (Cambridge  University  Press, 
Teddington,  England,  1952). 

E.  R.  S.  Winter,  Trans.  Faraday  Soc.  46,  81  (1950). 

(i2)  E Whalley  and  W.  G.  Schneider,  J.  Chem.  Phys.  20,  657  (1952). 

B.  N.  Srivistava  and  M.  P.  Madan,  Phil.  Mag.  43,  968  (1952);  Proc 
Phys.  Soc.  (London)  A66,  277  (1953). 


(14) 


Amdur,  Ross,  and  Mason,  J Chem.  Phys.  20,  1620  (1952). 
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suitably  described  by  an  exponential  form^^’^^ 
chosen  the  form 


Accordingly  we  have 


(p(r)* 


6 

I- Vet 


u <*0-r/rm) 

/ ¥'  \ k " 

- e 

1 -E2  I 

L* 

\ r ) 

(1) 


where  (p  (r)  is  the  potential  energy  of  the  molecules  at  separation  distance 

r , - 6 is  the  minimum  potential  energy,  rm  is  the  value  of  r for  which 

q>  (r)  is  a minimum,  and  K is  an  additional  parameter  which  may  be  considered 

a measure  of  the  steepness  of  the  repulsion  energy.  The  next  higher  term  in  the 

o 

dispersion  energy,  a term  varying  as  r and  representing  the  dipole-quadrupole 
energy,  has  not  been  included  since  the  same  effect  as  adding  such  a term  may 
be  achieved  to  a good  approximation  by  a variation  of  the  parameter  0 c . 


<15)  C.  Zener,  Phys.  Rev.  T7,  556  (1931). 

R.  H.  Fowler  and  E . A.  Guggenheim,  Statistical  Thermodynamics 
(Cambridge  University  Press,  Teddington,  England,  1939),  pp.  276-279, 
291-294,  give  a summary  of  the  evidence  prior  to  1939. 

R.  A.  Buckingham  and  J.  Corner,  Proc.  Roy.  Soc . (London)  A 189,  118 
(1947);  J-  Corner,  Trans  Faraday  Soc.  44,  914  (1948). 

<18)  I.  Amdur,  J.  Chem.  Phys.  _17,  844  (1949) 

(19*  J.  L.  Yntema  and  W.  G Schneider,  J.  Chem.  Phys.  18,  646  (1950). 

M.  Kunimune,  J.  Chem.  Phys.  18,  754  (1950);  Prog.  Theor.  Phys  5, 
412  (1950). 


(20) 
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A typical  potential  energy  curve  is  shown  in  Figure  1.  A more 

easily  visualized  measure  of  the  steepness  of  the  repulsion  energy  than  oc 

is  the  ratio  CJ  / r , where  0"  is  the  value  of  r for  which  CP  (r)  = 0 . 

Table  I gives  U /r  as  a function  of  at  ; included  for  comparison  is 

O’  / r for  the  Lennard  Jones  (12-6)  potential, 
m 

Table  I.  Ratio  of  the  low  velocity  collision  diameter,  (T  , to  the  position  of 

the  minimum,  rm  . 


Ot 

cr/r 

m 

12 

0.8761 

13 

0. 8832 

14 

0. 8891 

15 

0.8942 

Lennard- Jones  (12-6) 

0. 8909 

A potential  of  the  form  of  Eq.  (1)  was  derived  theoretically  by  Slater  and 

(21) 

Kirkwood'  ' for  the  special  case  of  helium,  but  its  use  as  a semi-empirical 
form  and  its  extensive  application  has  been  carried  out  mainly  by  Buckingham^ ^ . 


J.  C.  Slater,  Phys.  Rev.  32,  349  (1928);  J.  C.  Slater  and  J.  G.  Kirkwood, 
Phys.  Rev.  _37>  682  (1931). 

xt.  A.  Buckingham,  Proc.  Roy.  Soc.  (London)  A 168,  264  (1938); 

H.  S.  W.  Massey  and  R.  A.  Buckingham,  ibid.  A168,  378  (1938); 

A169,  205  (1938);  Buckingham,  Hamilton,  and  Massey,  ibid.  A 179, 

103  (1941);  R.  A.  Buckingham  and  R.  A.  Scriven,  Proc.  Phys.  Soc. 
(London)  65,  376  (1952). 


Fig.  1.  A typical  intermolecular  potential  energy  curve. 
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This  potential  has  the  defect  that  at  a small  separation  distance,  r , it 

c max 

spurious 

has  a/maximum,  and  approaches  minus  infinity  as  r approaches  zero  There- 
fore Buckingham  usually  has  used  Eq  (1)  only  for  r > r , and  for  r<  r 

m m 

used  another  form  in  which  the  inverse  sixth  power  term  of  Eq.  (1)  was 
multiplied  by- another  exponential  term  to  prevent  the  potential  from  having 
a maximum  However,  the  maximum  occurs  at  such  high  energies  that  it  has 
little  effect  on  ordinary  thermal  collisions,  and  we  believe  it  is  preferable 
to  use  Eq.  (1)  for  r > rmax  > and  to  set  the  potential  equal  to  infinity  for 
r<  rmax  . ■ We  have  called  this  form  a modified  Buckingham  potential,  which 
might  be  referred  to  as  the  "Exp-Six"  potential. 


H.  GENERAL  FORMULAS 

According  to  the  theory  of  Enskog  and  Chapman,  the  transport 
phenomena  depend  on  the  intermolecular  potential  through  collisions  of  single 
pairs  of  molecules,  and  the  only  feature  of  a collision  which  it  is  necessary  to 
consider  is  the  angle  of  deflection,  X , through  which  the  relative  velocity 
vector1  of  a pair  of  molecules  under  consideration  is  rotated  by  the  collision. 
The  angle  X is  given  by^^ 


X (v,b)  = 7T-  2b 


OO  _J_ 

[ I - b /r 2 - V(r">^uvz  ] dr/r2  , 


r 

c 


{!) 


E.  H.  Kennard,  Kinetic  Theory  of  Gases  (McGraw-Hill  Book  Go.,  Inc 
New  York,  1938),  pp  115-122. 

Reference  (9b),  Chapter  1. 
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■where  v is  the  absolute  value  o£  the  relative  initial  velocity  of  the  pair  of 
molecules,  b , the  collision  parameter,  is  the  perpendicular  distance  between 
one  of  the  molecules  and  the  initial  line  of  relative  approach  of  the  other,  u, 
is  the  reduced  mass  of  the  colliding  pair,  and  r£  is  the  distance  of  closest 
approach  during  the  collision,  given  by 

/-  b Vrt4-  2Cj)(rc)/Jay2=  o . 


(3) 


The  reduced  velocity- dependent  cross  sections  are  determined  from  X.  by 


Air  . , n-(-Dlr' * 

S (K)=ai-a7Ti  a-cosLx>/6<i/s  , 

* •'n 


(4) 


m 


where  we  have  introduced  the  reduced  quantities  Ky/.v  /2fc  and  p =b/r 
Finally,  the  reduced  collision  integrals,  which  directly  determine  the  transport 
properties,  are 


<!,r0»  * 

SI 


(n+-i)  •*  ( T *)' 


o<-a 


eK/T  K^'s^OOdK,  (5) 


where  T*=  kT/e  , k is  Boltzmann's  constant,  and  T is  the  absolute 
temperature.  The  cross  sections  have  been  reduced  by  dividing  by  Tf  i‘ 


U) 


m 


and  the  collision  integrals  are  just  the  SI  (n.)  of  Chapman  and  Cowling^®), 


divided  by  their  values  for  rigid  spheres  of  diameter  r 
and  ^ (T  ) are  unity  for  rigid  spheres. 


m 


Thus  S 


a ) 


(K) 


(25) 


Reference  (4),  pp.  157,  160. 
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The  values  of  the  transport  coefficients  may  be  expressed  as 
infinite  series,  and  higher  approximations  to  the  coefficients  are  obtained 
the  more  terms  of  the  series  that  are  taken.  Fortunately,  the  convergence 
is  rapid,  and  very  few  terms  are  needed.  The  coefficient  of  viscosity,  71 
of  a pure  gas  is 


(6) 


where  M is  the  molecular  weight,  R is  the  gas  constant  per  mole,  and 
f ^ represents  the  infinite  series,  and  is  a complicated  function  of  the 
collision  integrals.  The  term  f ^ is  nearly  unity,  and  is  a slowly  varying 
function  of  T*  ; in  Appendix  I the  expression  for  the  first  three  terms  of  f^ 
is  given  in  terms  of  the  collision  integrals. 

The  coefficient  of  thermal  conductivity,  A , of  a pure  gas  is 


fx 

n(a'a)*cn 


(7) 


where  Cy  is  the  molar  heat  capacity  at  constant  volume  of  the  gas,  and  f ^ 
is  a term  similar  to  f ^ . The  first  three  terms  of  f ^ are  given  in  terms 

of  the  collision  integrals  in  Appendix  1.  The  formula  (7)  does  not  take  account 
of  the  large  amount  of  energy  transported  by  the  internal  degrees  of  freedom 
of  the  molecules,  and  thus  applies  only  to  monatomic  gases.  An  approximate 
correction  for  this  extra  energy  transport  may  be  made  by  multiplying  the 


* 


\ 
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right  hand  side  of  Eq  (7)  by  the  "Eucken  factor i<(^)(27)(28)^  which  is  usually 
taken  to  be  (9  V / 10  - 1/2)  , where  Y is  the  ratio  of  the  constant  pressure  heat 
capacity  to  the  constant  volume  heat  capacity  of  the  gas  No  exact  method  for 
making  this  correction  has  yet  been  worked  out. 

The  coefficient  of  mutual  diffusion  of  a binary  mixture  of  gases  of 
types  1 and  2 is 

n . j_  ( M'  i M*  3rt  \a  fp , 

UI2  ft  \ M,Ma  TT  / n (rj^  (8) 

•where  Dj^  i®  the  mutual  diffusion  coefficient,  and  are  the  molecular 

weights  of  the  species  1 and  2,  and  n is  the  total  molecular  density.  The 
subscripts  1,  2 refer  to  an  interaction  between  a molecule  of  type  1 and  one 
of  type  2 . The  term  f^  is  a very  complicated  function  not  only  of  T*  , but 
also  of  hC  , M , and  the  composition  of  the  diffusing  mixture.  The  dependence 
of  on  the  composition  is  entirely  contained  in  this  term,  as  well  as  the 

dependence  of  D ^ on  1-1  and  2-2  molecular  interactions.  The  third 
approximation  to  f^  is  given  in  Appendix  I.  It  seemed  worthwhile  to  extend 
the  expressions  for  the  diffusion  and  thermal  diffusion  coefficients  to  one  more 

A.  Eucken,  Phys.  Zeit.  14,  324  (1913). 

(27)  Refereilce  (4)(  pp,  237-240. 

(28)  Reference  (9b),  Chapter  7. 
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approximation  than  has  previously  been  used  for  a realistic  intermolecular 

(29) 

potential,  inasmuch  as  there  is  evidence'  ' that  the  lower  approximations 
might  differ  from  the  limiting  theoretical  values  by  amounts  greater  than  the 
error  involved  in  determining  the  quantities  experimentally.  The  expression 
•.or  the  self- diffusion  coefficient,  , is  obtained  from  Eq.  (8)  by  setting 

M = and  deleting  the  subscripts  1,  2. 

The  general  expression  for  the  thermal  diffusion  ratio  is  quite 
complicated,  and  for  computational  purposes  is  perhaps  best  expressed  in  the 
determinant  notation  of  Chapman  and  Cowling.  In  this  form,  the  expression 
for  the  gi-th  approximation  to  the  thermal  diffusion  ratio,  [ kjl  m » of  a 
binary  mixture  is 

(m.)  M,  -*•  , <m>,  M . + M,  \^a 


M 


cm.),  n,  + n,y«a  , n + n2 

■4, 


f oi  \ •crjf  0/-1  ’ ^ 

TJm  p-  j$lnx) 


00 


(9) 


jf  (m) 

where  x^  and  x-j,  are  the  mole  fractions  of  the  two  components,  dfr  is 

a determinant  of  (2m  + 1)  order,  whose  general  term  is  a..  , where  i and 

/tim) 

j range  from  -m  to  m . is  the  minor  of  cv  obtained  by 

striking  out  the  row  and  col.  nn  containing  a_  . Th'i  elements  a^  are 
functions  of  the  collision  integrals  »n)*,  the  molecular  weights  of 

the  component  gases,  and  the  composition  of  the  mixture.  The  egressions 
for  the  a.,  for  m - 2 are  given  in  Appendix  I. 

* J 


(29) 


Reference  (4),  pp.  169,  196. 
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For  binary  mixtures  of  heavy  isotopes,  where  the  forcf  ” between 

the  different  molecular  species  are  the  same,  Eq.  (9)  may  be  put  into  a 

simpler  form  by  expanding  the  determinants  in  powers  of  the  small  quantity 

(M  - M ) /(M  + M ) , and  keeping  only  the  first  power.  To  this  order  of 
12  12 


approximation, 


(10) 


where  k*  I is  a reduced  thermal  diffusion  ratio,  and  is  a function  of  the 
reduced  temperature  alone.  Since  the  expression  for  [k*J  m in  terms  of 
the  collision  integrals  is  complicated,  it  is  placed  in  Appendix  I.  We  prefer 
the  function  k^,  to  the  more  frequently  encountered  , for  reasons  which 

are  discussed  in  detail  at  the  end  of  Appendix  I.  The  two  functions  are  very 
simply  related:  =(118/10  5)  k*  . 

For  mixtures  of  gases,  the  following  three  ratios  of  collision 
integrals  occur  frequently: 


A*  = 


B*  = 


ia  ' 12  5 

[5Xi,a  -4na  ]/o,. 


0,1)  * 


(11) 


^ o.a)*  / 0,0  * 

c*  = n /a 

i a ' <7. 


These  quantities  are  equal  to  unity  for  rigid  spheres. 
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HI.  NUMERICAL  CALCULATIONS 

The  first  step  in  the  evaluation  of  the  collision  integrals  is  the 
calculation  by  numerical  integration  of  the  angle  oi  deflection,  X , as  a 
function  of  v , b , and  the  potential  parameters  Eq.  (2)  is  inconvenient 
for  numerical  computations  because  the  integrand  becomes  infinite  at  the 
lower  limit  of  the  integral.  This  singularity  may  be  removed  by  substituting 
into  Eq.  (2)  the  value  of  b from  Eq.  (3),  and  making  the  change  of  variable 
sin  0 - r jr  . after  which  Eq  (2)  may  be  rearranged  into  the  form 

A/  2/3  ( s , sin  QtSin  e (Wa)e  f e sin  e - e 

TT  - i ~ : — l 3 

y J KO-C/tOy^  K(l~C,/&)  \ cos*e 


\ 
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The  numerical  integration  of  Eq.  (12)  was  sufficient  to  determine  X to  a few 

ten  .thousandths  of  a radian.  Values  of  X , together  with  the  corresponding 

values  of  /5  , have  been  tabulated  as  a function  of  y and  K for  & =12, 

13,  14,  and  15;  the  length  of  the  tabulations  prevented  their  publication 

in  the  present  paper^^.  The  tabulations  cover  thirteen  values  of  K for 

each  value  of  oc  , as  in  Table  III,  and  an  average  of  about  20  valueB  of  y for 

each  K .'  The  size  of  the  intervals  on  y varies  from  0.005  to  0.40,  depending 

on  the  behavior  of  X > and  was  chosen  so  as  to  obtain  an  accuracy  of  0.  1 

percent  in  the  integration  for  ^ ^(K)  . 

The  calculation  of  X was  complicated  by  the  fact  that  for  some 

collisions  which  occur  at  less  than  a certain  critical  energy,  K , the 

o 

molecules  may  "orbit l*  around  each  other  for  an  indefinite  number  of 
rev©  utions,  and  X consequently  tends  to  minus  infinity.  Such  orbiting 
collisions  have  been  discussed  in  connection  with  the  Lennard- Jones  (12-6) 
potential^^ ^ , and  the  same  general  behavior  is  observed  with  the  present 
potential.  For  energies  equal  to  or  less  than  Kq  , there  is  a critical  value 
of  beta,  /3  ^ , for  which  orbiting  will  occur.  There  are  two  values  of  the 
distance  of  closest  approach,  y , which  correspond  to  /3  ; values  of  y 

between  these  two  critical  ones  cannot  occur  physically.  The  larger  y , 
designated  as  y^  , will  cause  the  integrand  of  Eq.  (12)  to  become  infinite  at 
the  upper  limit,  so  that  the  value  of  yQ  may  be  determined  by  numerical 

(30v 

* ' Multig&phed  tables  may  be  obtained  from  the  author  on  request. 
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solution  of  the  equation 


{(d/oC  )g, 


& ( i - yo ) 


K(l-  (t/oc)  y 1 


K (\-W*) 


-n-o 


(15) 


which  is  obtained  from  Eqs.  (12)  and  (14).  The  value  of  A may 

' o 

be  found  by  substituting  the  value  of  y into  Eq.  (13)  and  solving.  The 

0 

smaller  value  of  y , designated  as  y , also  satisfies  Eq.  (13)  with  6 - P>  , 

1 • ' ' o 

and  may  be  calculated  in  this  way.  At  y = v^  the  integrand  of  Eq.  (12)  becomes 

infinite  for  (p  = arc  sinfy^/y^.  Between  y = y^  and  y = y^  there  are  two 

values  of  y which  give  the  same  when  substituted  into  Eq.  (13),  but 

these  values  are  without  physical  significance.  Values  of  * y , and  y^ 

for  several  values  of  K are  given  in  Table  II;  y and  y become  equal  at 

o 1 

K = K . 
o 

For  small  values  of  X > if  is  useful  to  have  the  asymptotic 
series,  valid  for  large  values  of  y : 


x-~- 


isn 

' 1 

51  S’  TT 

1 

1 (o 

. K(  i - <o/cl)  y6 

lOi* 

- 

(6/o t ) e 


o*(i-y) 


K ( I - 4/0' ) 


V 2 j Uy 


UfatK/)3 

v 


(16) 


which  was  obtained  by  the  method  of  Amdur  and  Pearlman 


(31) 


I.  Amdur  and  H.  Pearlman,  J.  Chem.  Phys.  % 503  (1941); 
I.  Amdur,  J.  Chem.  Phys.  15,  482  (1947). 


I 


Table  II.  Conditions  i 


K 

A 

12 

.4 

1.  7916 

K = . 78101 
o 

1. 5835 

13 

.4 

1. 7716 

. 8 

1. 5636 

K = .82691 
o 

1. 5530 

14 

.4 

1. 7555 

.8 

1. 5530 

K = . 86809 
o 

1.  5281 

15 

.4 

1.  7421 

.8 

1. 5437 

K =.90530 
o 

1.  5075 

! 

r 


WIS-ONR- 1 
19  June  1953 

16 


orbiting. 


Vo 

y, 

1.4468 

.9694 

1.  1849 

1. 1849 

1. 4853 

.9649 

oo 

O' 

fH 

1.  1076 

1.  1708 

1.  1708 

1. 4253 

.9621 

1. 2209 

1.0710 

1.  1594 

1. 1594 

1.4166 

.9604 

1. 2248 

1.0504 

1.  1498 

1. 1498 
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To  evaluate  the  cross  sections 


integration  variable  was  changed  from  j3 


I 5 ) 

Sv  (K)  defined  in  Eq  (4),  the 
to  y , and  Eq.  (4)  transformed  to 


(l) 

S («)■ 


2 

r 

(1-  COS'**.) 

' 1 l i + c-D1! 

* 1+1 

J 

k 

K(l-C/oOy4  K (I- 4/oc) 


tx(i-y) 

+ /(Ly-,') 

u f 1 - t /rJ\  \ : ' 


y°y 


(17) 


where  yc  i®  t^Le  value  of  y for  which  = 0 , and  represents  the  reduced 
distance  of  closest  approach  for  a head-on  collision.  Most  of  the  integration 
could  then  be  carried  out  numerically,  using  the  values  of  X tabulated  as  a 
function  of  y for  a given  K and  . To  avoid  the  necessity  of  calculating 
accurate  values  of  y^  from  Eq.  (13),  the  numerical  integration  was  started 
from  a lower  limit  y . which  was  slightly  greater  than  y , and  the  integration 

ai  C 

from  yc  to  y^  executed  as  follows.  Eqs.  (12)  and  (13)  were  differentiated 

to  form  (d  X /d / ) in  terms  of  X and  for  the  limit  of  — * 0 . 

Solution  of  the  resulting  differential  equation  then  showed  that  X varied  linearly 

with  when  ft  was  near  zero,  so  that  the  integration  from  y^  to  y^  could  be 

performed  analytically.  The  equations  obtained  for  these  small  contributions 

to  the  cross  sections  are: 

2. 


5ln(n-x‘)> 


(2)  2 2 3 /3. 

A (K)*  *f  + i (rT - xj 


- COS  l(  tT-XA) 


m. 

- -—A-  elf l 2 (tr-  X.)  , 
4 rr-XA 


(18) 


u ^ 

A0>(K)-/f-  1 F~Vy2.  fl-C0S3(jr-  X )]  - I /-A-  r,-C05(iT-  X 

'A  V"'V  L A - (tt  - X ) L * 

?TjiA  5ln  ! S7x  SinfK-XJ  , 

A A 
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A }(K)-  "I  + (^7-)3  [i-  cos  4(TT  - X,)] 


A 25$, 
2 


+ s -4_ 

16  (tt-X/ 


COS  2 (it- x 


'll  - — ' 
' J ~ tr- 


ain 4(tt-xJ 


(18) 
(cont.  ) 


5 & 

8 TT-Xa  u V ) 

where  the  ^ )(K)  are  the  contributions  to  the  ^(K)  for  the  range  yc 

to  yA  , and  /S  A and  X A are  the  values  of  /i  and  X corresponding  to 

y . The  numerical  integration  of  Eq.  (17)  when  no  orbiting  occurred 
A 

(K  > K ) was  carried  out  from  yA  to  y^,  , where  y^  was  a large  value 
of  y , and  the  remainder  of  the  integration,  from  y to  infinity,  carried 

2*4 

out  by  substituting  the  small  angle  formula,  Eq.  (16),  into  Eq.  (17),  expanding 
(1  - cos  ^ X ) » and  integrating  analytically  to  yield: 


If,  1 n-t-'l1 * 1 

FU)(K)  = 

2 

,a 

■ 1 ( « ^ 

^ 23  ( - \3] 

J.I'  a.  ,+i  J 

i 

• 

14  ) 

yE 

io\K0-i/K)y‘ J 

512  VK(i-4/*)y£  j . 

oc  O'V  1 

( !£H)  fc/oQe  J_ 

V U>)  K ( i - L/u) 


/ 2TT  \ar £ 

\ j L 8oc 


32focy): 


where  ^(K)  is  the  contribution  to  ^(K)  for  the  range  y^,  to 


infinity. 


i 

t 
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When  K is  equal  to  or  less  than  the  critical  energy  for  orbiting, 
Kq  , numerical  integration  of  Eq.  (17)  cannot  be  carried  out  over  the  entire 
range  from  y to  y because  (1)  the  region  from  y to  y is  physically 
inaccessible,  and  in  fact  X is  imaginary  there;  (2)  the  integrand  oscillates 
violently  near  y^  and  , where  X goes  to  minus  infinity.  For 
K < Kq  the  integral  from  y^  to  y was  therefore  broken  up  into  four 
integrals:  B<  i >(K)  , C(  1 >(K)  , D{  ^ *(K)  , and  E^  >(K)  . The  integral 
^ )(K)  extended  from  y^  to  y^  , where  y^  was  slightly  less  than  y^  , 

( l ) 

ahd  was  evaluated  numerically.  Similarly,  the  integral  E (K)  extended 

from  y^  , where  y^  was  slightly  greater  than  y , to  y^  , and  was  also 

evaluated  numerically.  The  integrals  C ^(K)  and  D^)(K),  cave  ring 

the  regions  y^  to  y^  and  yQ  to  y^  , respectively,  were  calculated  by  the 

(7) 

method  of  Hirschfelder , Bird,  and  Spotz'  ' , in  which  X is  represented 
near  the  regions  of  orbiting  by  equations  of  the  form 

X - cons-h/ (/*  V/)  , 


(20) 


and  Eq.  (4)  is  integrated  analytically  between  the  proper  limits.  The  results 


of  this  integration  are: 
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C°'(K)*(/V/361)[  Cos  x 0 + X B sin  XB  - X g CL  (-  Xg)]  , 

OO-  § (Aa-^  )L  I-  COS  2 x6  1-  2X6  S'O  nc  - t*;  C;  (-axj] 


20 


6 • > 


<8> 


C" «)• 


4-  cos  3 Xg  - 3 cos  Xg  + 3 Xe  sia-3X6 


+ 3X6slrX6-  r«„)-3X‘ci(-x,) 


(21) 


4),  5/2  i\ 

c (k)  - 2i(/50  -/*,  ) 


2 . 1\ 


5"  ~ COS  4 X - f COS  2Xa  +4.  X Sin  ‘h  X 
o o 0 o 

+ 8 XB  sin  X1.  c;  ('fXj 


-itX6  c.  (-ax0) 


where  and  X n are  the  values  of  and  X corresponding  to  yR  , and 


(32) 

Ci(-  X ) is  the  cosine  integral,  a tabulated  function'  , defined  as  folio 


ws 


Cl  (i) 


cos 


dt  . 


(21b) 


The  ^(K)  are  of  the  same  form  as  the  ^ ^(K). 

t > i \ 

The  cross  sections  S'-1-  '(K)  are  not  convenient  to  use  directly, 

since  they  vary  rapidly  when  K is  small,  and  tend  to  infinity  as  K goes  to 

(5) 

zero.  Kihara  and  Kotani  have  shown  that  in  the  limit  K — * 0 , the  cross 


Federal  Works  Agency,  Work  Projects  Administration  for  the  City  of 
New  York,  Tables  of  Sine,  Cosine,  and  Exponential  Integrals,  Vol.  II 
(1940);  Table  of  Sine  and  Cosine  Integrals  ( 1942) . 
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sections  are  independent  of  the  potential  energy  of  repulsion,  and  depend  on  the 

attraction  energy  alone.  From  this  it  is  easy  to  show  that  the  ^ ^(K)  are 

L / 3 

proportional  to  [k(1  - 6/c*  }]  for  K — * 0 , we  therefore  define  the 

functions  ^ ’(K)  : 


R(  1 ^K)  = [k(1-6/oi)]  1°S^)(K)t 


(22) 


which  are  given  in  Table  III  , and  plotted  for  ci  =12  and  15  in  Figure  2 for 
the  two  most  important  cases,  1 = 1 and  2.  For  comparison,  the  analogous 
functions  for  the  Lennard- Jones  (12-6)  potential  are  also  plotted.  Since 
R^  ^(O)  is  independent  of  the  repulsion,  the  values  given  in  the  table  for  1=1 
and  2 were  taken  directly  from  Kihara  and  Kotani^  . For  j - 3 and  4 , 
the  values  of  R'  (0)  were  obtained  by  extrapolation,  a procedure  which  is 
permissible  because  these  functions  occur  in  the  expressions  for  the  transport 
properties  only  as  correction  terms. 

For  use  in  evaluating  the  collision  integrals,  the  cross  section 
functions  were  approximated  by  the  following  algebraic  expressions; 

0 6 K < 0.2  • 

0.2  <K  < 2 , 


_Cl) 

(l)  / 

U) 

i 

, 3 

(t)  1 . 

II 
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(i)  v 

(t)  ly  a 
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a) 

a)..-' 

(i) 

5 (K)  = 

q-VK 

+ af 

K + 
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(0 

(0 
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co  -t 

0)  _4 

S 00- 
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& +■  0. 
9 

K 

<0 

+ an  K 

(jl) 
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-1 
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a«  * 
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«) 
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.1 
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-f 
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' + 0-(5  K 
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2<  K < 10  • 


10  6 K < °°  ■ 


(23) 
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Table  III  Cross  section  functions, 


(K)  - 


* 

- 

e 


R 


12 

13 

14 

15 

0 

1. 5778 

1. 5778 

1 5778 

1.5778 

0.4 

1 . 570  1 

1.  5588 

1 . 549  5 

1 . 5490 

0.78101 

0.8 

1. 5327 

1. 5280 

1.  5277 

| 

1. 5275 

1.0 

1.5178 

1.5297 

1. 5320 

1.5282  j 

1.  2 

1.4083 

1.4392 

1.4566 

1.4811 

1.6 

1. 2271 

1. 2714 

1.  30  27 

1. 3282 

2 

1. 1336 

1.  1743 

1 2073 

1. 2350 

3 

1.0361 

1.0805 

1. 1177 

1. 1509 

4 

1.0119 

1.0613 

1.  1018 

1. 1366 

5 

1 0109 

1.0626 

1.  10  70 

1. 1443 

10 

1.0646 

1.1297 

1.  1865 

1. 2356 

20 

1. 1546 

1. 2401 

1.3107 

1.3769 

50 

1.2824 

1.4032 

1.  5069 

1. 6000 

100 

1.3676 

1.5226 

1. 6579 

1.7811 

R*2)(K) 

0 

1. 7865 

1. 7865 

1 7865 

1. 7865 

0.4 

1.7923 

1.7936 

1. 7985 

1. 7980 

0.  78101 
0.  8 

1.7991 

1. 8001 

1.  3060 

1. 8045 

1.0 

1. 8357 

1 7864 

1.  7959 

1.7984 

1.  2 

v _ S941 

l on  ->n 

X . U IL  1 

1 n £ -i  —r 

X . OU4.  i 

1 A O A 1 

J.  . 

1.6 

1.8344 

1 8688 

1.  8826 

1.3918 

2 

1.6929 

1 7333 

1 7587 

1.7726 

3 

1 4442 

1 4921 

1. 5288 

1.5621 

4 

1 . 3418 

1. 3870  ' 

1 430  5 

1.4619 

5 

1.  2:986 

1.  3489 

1.  3918 

1.4274 

10 

1. 3051 

1. 3676 

1 4200 

1.4689 

20 

1.4141 

1.4972 

1.  5605 

1. 6272 

50 

1.6063 

1. 7233 

1. 8294 

1.9171 

100 

1. 7599 

1.9129 

2. 0461 

2. 1613 
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sections  are  independent  of  the  potential  energy  of  repulsion,  and  depend  on  the 

( n \ 

attraction  energy  alone.  From  this  it  is  easy  to  show  that  the  S ' '(K)  are 

1 j 3 

proportional  to  [k(1  - 6/ « ) ] for  K ~ * 0 ; we  therefore  define  the 

functions  R1^  ^ *(K)  r 

R'  £ *(K)  = [k(1  - 6/ot)J  (22) 


which  are  given  in  Table  IU  , and  plotted  for  OC  =12  and  15  in  Figure  2 for 
the  two  most  important  cases,  i = 1 and  2.  For  comparison,  the  analogous 
functions  for  the  Lennard- Jones  (12-6)  potential  are  also  plotted  Since 
r(  £ ){0)  is  independent  of  the  repulsion,  the  values  given  in  the  table  for  £ = 1 
and  2 were  taken  directly  from  Kihara  and  Kotani^  . For  i =3  and  4 , 

( i ) 

the  values  of  R (0)  were  obtained  by  extrapolation,  a procedure  which  is 
permissible  because  these  functions  occur  in  the  expressions  for  the  transport 
properties  only  as  correction  terms. 

For  use  in  evaluating  the  collision  integrals,  the  cross  section 
functions  were  approximated  by  the  following  algebraic  expressions: 

o -z  K £ o.8  j 
0.i<K  < 2 , 
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10  6:  K < °°  ' 
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Table  III.  Cross  section  functions. 


R(  ^ }(K) 


Continued 


r/3)(k) 


12 

13 

14 

15 

0 

(1. 8236) 

(1  8236) 

(1. 8236) 

(1. 8236) 

0.4 

1.8115 

1. 800  5 

1. 7938 

1. 7987 

0. 78101 

1. 7831 

0.8 

1 7749 

1 7821 

1. 7745 

1.0 

1.7796 

1. 7833 

1. 7778 

1. 7797 

1.2 

1. 6788 

1. 6930 

1. 7042 

1. 7285 

1.6 

1. 6121 

1. 6333 

1. 6433 

1. 6595 

2 

1. 5325 

1.5665 

1.  5884 

1. 6151 

3 

1. 3587 

1.4012 

1.4376 

1.4692 

4 

1. 2689 

1 3167 

1. 3557 

1.  3934 

5 

1 2256 

1. 2776 

1. 3206 

1. 3574 

10 

1.2222 

1. 2875 

1 . 3447 

1. 3893 

20 

1. 3151 

1.  3993 

1.4632 

1. 5363 

50 

1.4922 

1. 6047 

1 7045 

1. 7928 

100 

1. 6115 

1.7641 

J 

1. 8980 

L_ 

2.0141 

R(4)(K) 

0 

(1.9100) 

(1.9100) 

(1  9100) 

(1.9100) 

0.4 

1.9335 

1.9208 

1.9225 

1. 9274 

0.78101 

1.9345 

0.8 

1.9178 

1 9221 

1.9586 

1.0 

1.9330 

1 9191 

1.9065 

1. 9107 

1.2 

2.0278 

1 9948 

1.9907 

1.9674 

1.6 

1.9600 

1 9878 

1. 9770 

1. 9983 

2 

1.8895 

1.9176 

1 9300 

1 9258 

3 

1. 6584 

1. 70  26 

1. 7363 

1. 7663 

4 

1.  5197 

1.  5623 

1. 6045 

1 6348 

5 

1.4433 

1. 4941 

1 5328 

1 £705 

10 

1. 3868 

1. 4514 

1 4987 

1. 5488 

20 

1.4838 

1.  5670 

1. 6252 

1. 6893 

50 

1. 6931 

1 8089 

1.9105 

1.9954 

100 

1. 8636 

2 0136 

2 1450 

2.  2558 
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This  curve-fitting  procedure  effectively  permits  the  extrapolation  of  the 

( Jt  ) 

S'  (K)  to  infinite  K , which  is  necessary  because  the  potential  energy 
function  Cp  (r)  of  Eq.  (1)  misbehaves  for  very  small  values  of  r by 
becoming  negative. 

The  final  integrations  with  respect  to  K to  obtain  the  collision 

integrals  defined  in  Eq.  (5)  could  be  carried  out  by  replacing  S' x ;(K)  in 

Eq.  (5).  by  the  algebraic  expressions  obtained  from  Eqs.  (22)  and  (23) 

For  K >0.8,  the  resulting  integrals  could  be  easily  evaluated  in  terms  of 

elementary  functions  and  error  integrals,  which  are  tabulated  functions.  For 

K < 0.8  , the  integrations  were  performed  numerically  for  large  T*,  where 

this  region  did  not  contribute  greatly  to  the  result,  and  for  small  T*  by 

(33) 

means  of  incomplete  gamma  functions  tabulated  by  Kotani  „ In  practice  it 
turned  out  to  be  easier  also  to  integrate  numerically  for  K > 0 8 when  T* 
was  large.  For  numerical  integrations,  i algebraic  expressions  in  Eq  (23) 

( St  ) 

were  used  only  as  interpolation  functions  to  obtain  values  of  S '(K)  at 
small  intervals  of  K . 

The  collision  integrals  are  not  especially  well  suited  to  tabulation 
because  their  rapid  variation  for  small  T*  makes  interpolation  difficult 
unless  the  tabulations  are  made  for  small  intervals  of  T*  A more  suitable 
function  for  tabulation  is 

1 n)(T*}  * [ T*(l  - 6/<*  )]  1/3  A U 3 n)(T*)  (2.4) 


(33) 


M..  Kotani,  Ptoc  Phys  -Math  Soc.  Japan  24,  76  (1942). 
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This  function  varies  relatively  slowly  with  T*,  and  unlike  ^ ' n)*(T») 

is  finite  auid  independent  of  the  parameter  oc  when  T*  = 0 . The  functions 

2 ’}  ^ * n^{T*)  are  tabulated  in  Table  IV  for  T*  from  0 to  200  . The  inherent 

computational  error  is  estimated  to  be  of  the  order  of  0.  2 percent,  except 

for  very  small  T*  , where  it  may  be  as  large  as  one  percent.  To  maintain 

internal  consistency,  more  figures  ate  given  in  Table  IV  than  is  strictly 

justified  by  the  estimated  error.  The  reason  for  the  lower  estimated  accuracy 

( 1 ) 

at  small  T*  is  that  the  cross  section  functions  R'  '(K)  seem  to  have  an 
oscillatory  behavior  when  K is  small,  but  no  attempt  was  made  to  follow  this 
behavior  in  detail  since  many  calculations  at  closely  spaced  values  of  K would 
have  been  necessary.  The  labor  involved  did  not  seem  worthwhile,  particularly 
since  the  neglect  of  quantum  corrections  at  low  K is  probably  a much  more 
serious  error. 

IV.  ADDITIONAL  TABULATIONS 

In  Table  V are  presented  the  functions,  f , which  give  the  higher 
approximations  to  the  transport  coefficients  for  pure  gases;  the  superscripts 
indicate  that  these  are  the  third  approximations  rather  than  the  exact  values. 

The  function  f for  a pure  gas  refers  to  the  coefficient  cf  self-diffusion. 

The  second  approximation  to  the  reduced  thermal  diffusion  ratio  for 
a binary  mixture  of  heavy  isotopes,  [kip]  ^ , is  given  in  Table  VI  as  a 
function  of  T*  , and  plotted  in  Figure  3 for  a =12  and  15.  For  comparison, 
the  function  [k*j^  for  the  Lennard- Jones  (12-6)  potential,  as  calculated  by 
the  method  of  Chapman  and  Cowling  is  also  shown  in  Figure  3.  It  is  of  interest 
to  note  that  [ k*  ] ^ is  always  positive  for  c.:  - 15,  showing  that  there  is  no 
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inversion  temperature  when  the  repulsive  portion  of  the  potential  (p  (r)  is 
steep  enough. 

Table  VII  gives  the  values  of  the  ratios  A*  , B*  , and  C* 
which  occur  in  the  expressions  for  the  transport  properties  of  mixtures 
The  dimensionless  group  ? D^/lj  ‘ ^re<lueat^y  encountered  in  kinetic 
theory,  is  given  by  (bt^/51^  )A* 

Experimental  viscosity  results  are  often  expressed  by  an  equation 

of  the  form 

| = cT8’.,  (25) 

where  c and  s ^ are  supposed  to  be  constants  for  a particular  gas.  The 
temperature  index,  « , is  seen  to  be  given  by 


s^  = d log  7|  /d  logT . (26) 

The  temperature  index  has  the  following  physical  significance:  if  the 
intermolecular  potential  were  of  the  form  (p  (r)  = const,  /r  v , then  s ^ 
would  be  constant  and  equal  to  | 1/2  + 2 /(  V - 1)^  . We  can  investigate 


the 


behavior  of  s^  for  the  potential  given  by  Eq.  (1)  by  requiring  that  the  values 

of  T)  and  d 7]  /dT  determined  from  Eq  (25)  agree  with  those  from  Eq  (6) 
l i 

at  a given  temperature.  In  this  way  we  find  that  s.  is  not  a constant,  but 


n 


the  following  function  of  T* 


9 

Y ^ 


4X1 


(3,3) 


.a 


( 3,3)* 


d (09  fq 


d I , 


(27) 


where  we  have  made  use  of  the  recursion  relation 


d XI 


fin)* 


Xn+  2) 


, (A.o-fO* 

a 


IOQ 


SI 


(l.n)  * 


- I 
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(28) 


Another  common  method  of  expressing  experimental  viscosity  results  is  by  the 
Sutherland  equation, 


c'  T a 

V / + s^/r 


(29) 


where  c'  and  are  supposed  to  be  constants  for  a particular  gas  This 

equation  is  based  on  the  physical  model  of  rigid  spheres  which  have  a weak 
attraction  for  one  another.  The  physical  significance  of  the  Sutherland 
constant,  S ^ , is  that  it  is  proportional  to  the  potential  energy  of  two  such 

spheres  when  in  contact.  We  can  find  the  dependence  of  S,  on  T*  in  the 


same  way  that  the  behavior  of  8..  was  found;  it  is 


Sn  _\ 

T 


XI 


0,3)* 


XI 


7517* 


1 

d /oq  T* 


!-  *(/- 


XI 


(3,3)* 


n<a*a>*  ) d 


dlcxj 


(30) 


An  entirely  analogous  result  holds  for  the  coefficient  of  thermal 
conductivity,  A - The  expressions  for  s and  S ^ /T  are  the  same  as 
given  in  Eqs.  (27)  and  (30),  respectively,  except  that  f ^ is  replaced  by  f ^ 
In  a similar  way,  we  can  define  a temperature  index  for  self- 

diffusion, 


®D  = d D11  /d  1o8  t » 


(31) 
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and  a Sutherland  constant  tor  self  - diffusion. 

U rp  7. 


D 


11 


VT 


(32) 


If  Dji  in  Eq.  (31)  is  expressed  at  constant  density,  we  find  that  the  temperature 
index  s^  is  the  following  function  of  T* 


s 

D 


1 

2 


3X1  _ cllo^iD 

SI0'"*  + d !o j T *■ 


(33) 


For  the  Sutherland  constant 


VT 


we  obtain 


a 


0,1)* 


a 


(1,1)4 


A'iikJ 

d !o(j  T*  J 


Ci,a>* 

Ci,  i ? * 


* loc) 

d locj  T * 


(34) 


independent  of  whether  D is  expressed  at  constant  density  or  constant 
pressure. 

Values  of  the  temperature  indexes  for  viscosity,  thermal 
conductivity,  and  self-diffusion  were  calculated  from  Eqs  (27)  and  (33) 
neglecting  the  terms  involving  f , and  are  given  as  a function  of  T*  in 
Tabic  VTJT A . Tbc  error  introduced  by  nsglscting  tbs  terms  in  f is  in  3.11 
cases  less  than  one  percent.  In  Figure  4A  is  plotted  sv  vs  T*  for 
°(  = 12  and  15  , together  with  the  corresponding  s ^ for  the  Lennard- 
Jones  (12-6)  potential  for  comparison.  Values  of  the  Sutherland  constants 
are  given  in  Table  VUIB,  and  plotted  in  Figure  4B  for  a = 12  and  15  , 
and  the  Lennard- Jones  (12-6)  potential.  The  neglect  of  the  terms  involving 
f in  the  calculations  of  the  Sutherland  constants  introduced  an  error  of  less 


than  one  percent. 
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Numerical  values  of  the  potential  parameters  of  Eq,  (1)  for  a 
particular  gas  may  be  calculated  by  comparison  of  the  experimentally 
determined  values  of  S/T  or  s with  the  theoretical  values  in  Table  Yin. 
Caution  must  be  exercised  if  this  procedure  is  used,  however,  because 
inconsistent  results  may  be  obtained  unless  the  experimental  quantities  are 
calculated  in  the  same  way  as  the  theoretical  quantities  were;  that  is,  from  both 
the  temperature  derivative  and  the  numerical  value  of  the  transport  coefficient 
at  a given  temperature. 
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Table  IV  Collision  IntegT  al  Functions  ^ 

za'")(T'') -[T  " (I-  T'-  JkT/e 


oc=  12 


zO.  i) 

zO,  2) 

zi1-  3> 

Z1-1  4> 

5>' 

Zt3>  3> 

0 

1.  1870 

1.0551 

0 9672 

0 9027 

0.8572 

1 1178 

. 1 

1.  1911 

1.0  523 

0.9584 

0.8883 

0.8320 

1.  1079 

. 2 

1.  1662 

1.0142 

0 9025 

0 8113 

0 7338 

1 0704 

. 3 

i .1243 

0.9515 

0 8232 

0 7239 

0.6475 

1.0161 

4 

1.0750 

0.8900 

0 7597 

0 6667 

0.6005 

0 9641 

, 5 

1.0  282 

0.840  2 

0,7160 

0.6328 

0.5764 

0.9194 

. 6 

0.9873 

0.8025 

0.6870 

0,6131 

0, 5643 

0 8828 

. 7 

0.9530 

0.7745 

0.6679 

0.6019 

0. 5589 

0.8536 

.8 

- 0.9248 

0.7538 

0.6555 

0.5959 

0.5572 

0. 8305 

.9 

0.9016 

0.7387 

0.6475 

0. 5930 

0.5576 

0. 8125 

1.0 

0.8825 

0.7272 

0.6425 

0. 5923 

0. 5593 

0.7985 

1.  2 

0.8541 

0.7130 

0 6385 

0 5942 

0,5646 

0 7793 

1.4 

0.8350 

0. 70  58 

0 6388 

0 5985 

0.5711 

0 . 7684 

1.  6 

0.8221 

0 7028 

0,6414 

0.6039 

0.5779 

0.7625 

1 8 

0.8135 

0 7025 

0.6451 

0.6096 

0 5845 

Or.  7601 

2.0 

0 8080 

0. 7037 

0.6496 

0.6155 

0. 5910 

0.7600 

2.  5 

0.8023 

0. 7104 

0.6616 

0. 6296 

0.6059 

0. 7654 

3.0 

0.8031 

0 7195 

0.6735 

0.6425 

0.6180 

0.7746 

3.  5 

0.8070 

0. 7288 

0 6847 

0,6541 

0.6305 

0.7853 

4 

0.8126 

0.7382 

0.6950 

0.6644 

0. 6404 

0.7963 

5 

0. 8253 

0.7556 

0 . 7129 

0.6819 

0.6572 

0.8174 

6 

0 8384 

0. 7709 

0.7280 

0.6962 

0. 6705 

0.8364 

7 

0 8509 

0.7842 

0.7409 

0. 7081 

0.6817 

0.8533 

8 

0.8623 

0.7961 

0 . 7520 

0 7ib5 

0.6912 

0 8680 

9 

0.8730 

0 8067 

0 , 7618 

0 7272 

0 6993 

0 8813 

10 

0. 8829 

0 8160 

0. 770  2 

0 7351 

0.7066 

0 8929 

12 

0.9002 

0 8322 

0.7848 

0.7485 

0 7189 

0.9133 

14 

0.9150 

0.8457 

0.7969 

0 7594 

0 7292 

0.9301 

16 

0 9280 

0.8572 

0 8074 

0. 7690 

0.7380 

0.9444 

18 

0 9394 

0 8674 

0 8164 

0. 7773 

0, 7459 

0.9568 

20 

0 9497 

0 . 8762 

0 . 8245 

0 . 7846 

0. 7530 

0 9680 

25 

0,9710 

0 8951 

0 . 8415 

0 8007 

0. 7682 

0.9910 

30 

0.9885 

0. 9103 

0.8555 

0 .8138 

0 7808 

1.0094 

35 

1.0032 

0.9232 

0 8674 

0 8253 

0 7921 

1 0250 

40 

1.0157 

0.9346 

0 . 8781 

0 8355 

0.8021 

1.0385 

45 

1 . 0 270 

0.9446 

0 . 8876 

0 8450 

0 8114 

1 0508 

50 

1.0371 

0.9538 

0 8965 

0,8535 

0 8199 

1.0614 

60 

1 . f 549  ** 

0.9701 

0.9123 

0,8691 

0. 8355 

1 0807 

70 

1.0700 

0.9846 

0.9264 

0.8827 

0 8495 

1 0978 

80 

1.0834 

0.9973 

0.9388 

0 8957 

0 8622 

1. 1129 

90 

1 0959 

1 0091 

0.950  4 

0.9074 

0 8743 

1 1265 

100 

1 1071 

1 0201 

0.9615 

0 9184 

0 8853 

1.  1391 

. 200 

1. 1906 

1 , 1038 

1.0462 

1 0049 

0 9729 

1. 2356 
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Table  IV  Collision  Integra.!  Functions  continued 


oC  = 12 


T 

z<2  2) 

.2.  4> 

Z 

'2.  5) 

z{2.  6) 

4) 

0 

1.  1947 

1.0951 

1 

0221 

0 

9706 

0 9193 

1.0927 

1 

1 1983 

1 1002 

1 

.0291 

0 

9 750 

0 ,9323 

1. 1062 

2 

1 2098 

1.  1094 

1 

.040  1 

0 

9 770 

0 9188 

1.1119 

3 

i 2085 

1.0936 

0 

9978 

0 

9129 

0.8381 

1.0843 

1 1799 

1 0471 

0 

9367 

0 

8453 

0 7711 

1.0  360 

5 

i, 1421 

0.9974 

0 

. 8835 

0 

7953 

0 7280 

0.9870 

.6 

1 1028 

0.9538 

0 

8428 

0 

7614 

0 . 70 1 7 

0.9450 

7 

1„8667 

0 9185 

0 

.8133 

0 

7391 

0.6862 

0 9111 

. 8 

1,0352 

0 890  7 

0 

7923 

0 

7249 

0 6775 

0 8847 

9 

1 0084 

0 8692 

0 

7775 

0 

7160 

0.6731 

0.8642 

1 0 

0.9859 

0 . 8528 

0 

, 7673 

0 

7108 

0.6716 

0.8488 

1.  2 

0,9520 

0 8309 

0 

7561 

0 

7075 

0 6736 

0.8282 

1.  4 

0.9290 

0 ,8189 

0 

7525 

0 

. 7093 

0 6790 

0.8170 

I 6 

0 9137 

0 8 130 

0 

7 531 

0 

7139 

0 6860 

0. 8120 

1 8 

0.9039 

0 8111 

3 

7560 

0 

7198 

0 6936 

0. 8107 

2 0 

0.8978 

0 8116 

0 

7604 

0 

7263 

0 7014 

0.8119 

2.  5 

0.8930 

0.8189 

0 

7742 

0 

7435 

0. 7202 

0.8207 

3.0 

0 8963 

0 8300 

0 

7891 

0 

7600 

0 7374 

0.8331 

3,5 

0 9033 

0 8423 

0 

80  35 

0 

7752 

0. 7529 

0.8466 

4 

0 9122 

0 8546 

0 

8171 

0 

7892 

0 7668 

0.8599 

5 

0.9313 

0 8778 

0 

8415 

0 

8135 

0. 7906 

0 . 8849 

6 

0,9503 

0 8987 

0 

8623 

0 

8339 

0. 8105 

0.9072 

7 

0.9681 

0 9172 

0 

8806 

0 

8516 

0.8275 

0.9269 

8 

0.9845 

0 9339 

0 

8967 

0 

8670 

0 8424 

0 9445 

9 

0 9995 

0 9488 

0 

9110 

0 

8808 

0 8559 

0 9602 

10 

1.0133 

0 9624 

0 

9241 

0 

8935 

0 8678 

0.9745 

12 

1 0381 

0 9862 

0 

9469 

0 

9153 

0 8891 

0 9998 

14 

) 0598 

1.0068 

0 

9664 

0 

9341 

0 9073 

1.0213 

16 

1..0788 

1 0 248 

0 

9838 

0 

9 508 

0.9236 

1 0404 

18 

1.0958 

1 0409 

0 

9993 

0 

9658 

0 9383 

1 0575 

■>n 

i n n 

i n cc7 

1 

n i o a 

A 

A -7  A / 

* * * * 

A W V ^ * 

\J 

7I7U 

U 7 3AU 

X . U (Ji 

25 

1 1448 

1 0875 

1 

0441 

1 

0095 

0.9812 

1.  1070 

30 

1 1729 

1 1145 

1 

0 70  3 

1 

0353 

1,0067 

1 1357 

35 

1 4976 

1 1382 

1 

0935 

1 

0582 

1 0 294 

1 1608 

40 

1.  2193 

1 . 1 59  3 

1 

1143 

1 

0787 

1 0499 

1. 1832 

*5 

1 2393 

1 1786 

1 

1332 

1 

0976 

1.0688 

1 2041 

^0 

1.2573 

1 1962 

1 

1509 

1 

1152 

1 0863 

1. 2228 

60 

1.  2901 

1 2283 

1 

1826 

1. 

1469 

1 1180 

1 2572 

70 

1 3189 

1 2571 

1 

2113 

1 

1756 

1 1468 

1. 2878 

80 

1 3451 

1 2828 

1 

2370 

1 

2014 

1 1730 

1. 3157 

90 

1 3690 

1 3068 

1 

2609 

1 

2257 

1 1973 

1 3413 

100 

1 3915 

1.  3292 

1 

2835 

1 

2481 

1 2198 

1.  3649 

200 

1 5600 

l 4988 

1 

4547 

1 

4 20  8 

1 3939 

1 . 5470 

1 
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i 


Ot*  13 


T* 

zU.  i) 

2) 

Z<1.  3) 

4) 

5) 

ZP.  3) 

0 

1. 1870 

1.0551 

0 9672 

0.9027 

0.8572 

1. 1178 

. 1 

1. 1813 

1.0447 

0.9528 

0 8846 

0 830  5 

1 1020 

. 2 

1 1618 

1.0150 

0 9084 

0.8217 

0.7479 

1 0718 

. 3 

1 1275 

0.9619 

0 8389 

0 7429 

0.6681 

1.0250 

. 4 

1.0851 

0.9071 

0.7807 

0.6893 

0.6236 

0.9787 

. 3 

1.0436 

0.8617 

0. 7400 

0 6575 

0 . 60 1 1 

0.9383 

.6 

1.0068 

0.8268 

0.7129 

0.6393 

0.5904 

0.9048 

. 7 

0 9756 

0.8009 

0.6953 

0.6294 

0 5862 

0.8780 

.8 

0.9498 

0 7818 

0.6841 

0.6244 

0. 5855 

0.8569 

.9 

0.9285 

0.7678 

0.6772 

0.6226 

0.5871 

0.8405 

1.0 

0.9112 

0. 7577 

0.6732 

0.6229 

0.5897 

0.8279 

1.  2 

0.8854 

0 , 7454 

0.6710 

0.6266 

0.  5970 

0.8111 

1.4 

0 . 8684 

0 . 7400 

0.6729 

0 6326 

0.6052 

0.8022 

1.6 

0.8573 

0.7386 

0.6770 

0.6395 

0.6136 

0.7981 

1.  8 

0.8503 

0. 7396 

0.6823 

0.6468 

0. 6219 

0.7974 

2.0 

0.8461 

0.7421 

0.6881 

0.6542 

0. 6299 

0.7987 

2.  5 

0.8436 

0.7519 

0.7032 

0.6715 

0.6483 

0.8073 

3.0 

0.8472 

0 7637 

0.7181 

0 6875 

0.6643 

0.8193 

3.  5 

0.8535 

0.7757 

0.7320 

0.7017 

0.6785 

0.8323 

4 

0.8614 

0.7874 

0.7446 

0 7145 

0. 6910 

0 . 8454 

■5 

0.8782 

0.8091 

0.7671 

0.7365 

0.7121 

0.8700 

6 

0.8950 

0.8282 

0.7862 

0.7547 

0. 7296 

0.8922 

7 

0.9107 

0.8451 

0.8024 

0.7702 

0 7441 

0.9118 

S 

0.9253 

0. 8601 

0 8166 

0.7836 

0.7568 

0.9294 

9 

0.9389 

0.8734 

0.8292 

0 7954 

0.7678 

0.9453 

10 

0.9514 

0.8855 

0. 8404 

0. 80  59 

0.7777 

0.9596 

12 

0.9736 

0.9065 

0. 8599 

0 8240 

0 . 7948 

0.9847 

14 

0.9927 

0.9245 

0.8765 

0.8394 

0 8094 

1.0063 

16 

1.0096 

0.9399 

0 8907 

0.8528 

0.8223 

1.0250 

18 

1.0247 

0.9537 

0.9034 

n OAAQ 

V • w*  u 

n on*? 

V*  , UJJ  1 

1 A 4 i n 

i . \j*±l  t 

20 

1.0384 

0.9659 

0.9147 

0.8754 

0 8440 

1 0569 

*5 

1.0674 

0.9922 

0.9392 

0.8987 

0.8664 

1.0895 

30 

1.0913 

1.0139 

0.9596 

0.9184 

0 8855 

1 1166 

15 

1. 1118 

1.0  328 

0. 9775 

0 9357 

0.9027 

1.  1403 

40 

1. 1296 

1.0492 

0.9933 

0.9510 

0.9179 

1.  1616 

45 

1. 1456 

1.0643 

1.0076 

0.9651 

0.9318 

1. 1806 

H 

1.  1602 

1.0778 

1.0208 

0.9783 

0.9447 

1. 1983 

60 

1.1861 

1.  1023 

1.0447 

1.0017 

0.9682 

1.  2297 

70 

1. 2087 

1. 1239 

1.0659 

1.0230 

0.9894 

1. 2577 

80 

1.2287 

1 1431 

1.0849 

1.0422 

1.0089 

1. 2830 

90 

1 . 2469 

1.  1608 

1.  1025 

1.0598 

1. 0267 

1. 3063 

100 

1.2639 

1 1774 

1.  1193 

1.0762 

1. 0434 

1. 3277 

200 

1.  3892 

1 3027 

1. 2456 

1. 2042 

1 . 1 723 

, 1.4915 

Table  IV. 


Collision  Integral  Functione  continued 


w*.  c:.:r-i 

19  June  1953 


a = 13 


T* 

ESI 

Z(2.  3) 

Z (2.  4) 

Z(2,  5) 

z(2,  6) 

Z<4>  4) 

0 

1.1947 

1.0951 

1 0221 

0 9706 

0.9193 

1.0927 

. 1 

1.  1985 

1.0994 

1 0269 

0.9709 

0.9265 

1.0986 

.2 

1.2052 

1 1090 

1 0376 

0.9783 

0.9247 

1. 1116 

. 3 

1. 2056 

1 0977 

1.0077 

0.9278 

0.8562 

1 . 0944 

4 

1. 1862 

1.0602 

0.9549 

0 8664 

0. 7934 

1.0526 

„ 5 

1. 1549 

1.0164 

0.9059 

0.8190 

0. 7519 

1 0076 

.6 

1.1207 

0.9766 

0.8676 

0.7866 

0.7267 

0 9681 

.7 

1.0882 

0.9436 

0.8395 

*0.7653 

0. 7120 

0.9360 

.8 

1.0594 

0.9175 

0.8196 

0.7519 

0 7043 

0.9110 

.9 

1.0346 

0.8972 

0.8057 

0.7439 

0. 7009 

0.8917 

fVI 

1.0138 

0.8820 

0. 7965 

0. 7398 

0.  7005 

0.8775 

0.9824 

0.8618 

0. 7870 

0 7382 

0 . 7044 

0.8590 

0.9613 

0.8514 

0. 7850 

0.7418 

0. 7118 

0.8499 

0.9475 

0.8471 

0.7871 

0.7481 

0.7206 

0.8466 

1.8 

0.9391 

0 . 8465 

0.7916 

0.7557 

0 , 7299 

0.  8469 

0.9344 

0.8484 

0.7976 

0.7638 

0.7392 

0.8496 

0.9327 

0.8590 

0.8148 

0 . 7844 

0.7615 

0.8617 

0.9388 

0.8731 

0.8326 

0.8040 

0. 7817 

0.8771 

0 . 9485 

0.8880 

0.8497 

0.8219 

0. 7998 

0.8931 

0.9597 

0.9027 

0 8658 

0.8381 

0.8161 

0.9087 

0.9831 

0.930  3 

0.8944 

0.8668 

0.8443 

0.9378 

1.0058 

0.9548 

0.9191 

0.8911 

0.8679 

0.9637 

1.0270 

0.9769 

0.9408 

0.9121 

0.8885 

0.9867 

8 

1.0464 

0.9964 

0 9600 

0.9309 

0.9066 

1.0072 

1.0643 

1-0143 

0.9773 

0.9477 

0.9230 

1.0257 

10 

1.0809 

1.0306 

0.9929 

0.9628 

0.9377 

1.0425 

12 

1. 1105 

1.0592 

1 . 0 20  7 

0.9898 

0.9641 

1.0725 

14 

1. 1363 

1.0843 

1.0449 

1.0133 

0.9871 

1 0984 

16 

1.  1595 

1.  1064 

1.0662 

1.0342 

1.0076 

1. 1216 

18 

1 l on  i 
* • *uv * 

1 1 O L A 

A • lt.U7 

i a a r « 

k . UQ3  i 

1.0530 

1. 0262 

1.  1422 

20 

1. 1991 

1. 1448 

1.  1033 

1.0704 

1.0432 

1. 1611 

25 

1.2406 

1,  1847 

1.  .1423 

1 . 1085 

1.0812 

1. 20  27 

30 

1.  2758 

1.2187 

1. 1757 

1. 1416 

1. 1138 

1.2382 

55 

1.3066 

1.2486 

1 2052 

1.  1709 

I. 1430 

1. 2694 

40 

1.5344 

1.2757 

1. 2320 

1.  1975 

1. 1694 

1. 2974 

45 

1.3595 

1.300  4 

1. 2565 

1. 2217 

1. 1937 

1 . 3230 

50 

1.3826 

1.3232 

1.2792 

1 . 2444 

1. 2165 

1 . 3466 

60 

1.4243 

1.3644 

1 3202 

1 . 2854 

1. 2577 

1 3896 

70 

1.4615 

1.4015 

1. 3569 

1 . 3224 

1. 2945 

1 . 4280 

80 

1.4951 

1 . 4348 

1.  3906 

1. 3559 

1.  3282 

1.4628 

90 

1.5262 

1.4656 

1.4215 

1 . 3S69 

1.  3595 

1 . 4948 

100 

1.5546 

1.4942 

1.4500 

1.4161 

1.  38S5 

1.  5248 

300 

1.7718 

1. 7128 

1 6699 

1 6376 

1. 6114 

1 7518 
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0L=  14 


— 

T* 

Z<1.  1) 

Z<!.  2) 

z(l.  3) 

z(1*  4> 

Z*1’  5) 

z(3,  3) 

. _ _ 

0 

1. 1870 

1.0551 

0.9672 

0.9027 

0.8572 

1. 1178 

. 1 

1. 1727 

1.0382 

0.9481 

0.8814 

0.8287 

1.0987 

. 2 

1. 1575 

1.0146 

0.9115 

0. 8281 

0.7571 

1,0729 

. 3 

1. 1286 

0.9684 

0 8497 

0.7566 

0.6836 

1.0308 

,4 

1.0914 

0.9191 

0.7963 

0. 70  69 

0. 6420 

0.9891 

.5 

1.0542 

0.8776 

0.7585 

0 6773 

0.6214 

0.9524 

.6 

1.0208 

0.8455 

0.7336 

0.6608 

0.6121 

0.9218 

.7 

0.9924 

0.8217 

0.7176 

0. 6521 

0. 6091 

0.8974 

.8 

0.9688 

0 . 8042 

0. 70  76 

0.6483 

0.6096 

0.878.1 

.9 

0.9495 

0.7917 

0.7018 

0 . 6476 

0.6121 

0.8633 

1.0 

0.9337 

0.7827 

0. 6989 

0. 6488 

0.6157 

0. 8519 

1.  2 

0.9106 

0.7724 

0. 6984 

0.6542 

0 . 6246 

0.8373 

4.4 

0.8956 

0.7686 

0. 7019 

0.6617 

0.6342 

0.8300 

1.6 

0.8863 

0.7687 

0.7074 

0 . 6700 

0. 6440 

0.8273 

1.8 

0.8808 

0.7710 

0.7139 

0.6784 

0 . 6534 

0.8277 

2.0 

0.8780 

0.7747 

0.7208 

0.6868 

0.6625 

0.8300 

2.5 

0.8784 

0.7871 

0.7385 

0.  7068 

0.6836 

0.8406 

3.0 

0. 8843 

0.8011 

0.7557 

0. 7252 

0.7022 

0. 8545 

3.5 

0.8928 

0.8152 

0. 7716 

0. 7416 

0.7188 

0.8692 

4 

0.9025 

0.8288 

0.7863 

0.7565 

0.7335 

0 . 8840 

5 

0.9228 

0.8539 

0 8125 

0. 7824 

0.7587 

0.9121 

6 

0.9426 

0.8764 

0 8349 

0. 8042 

0. 7796 

0.9376 

7 

0.9612 

0.8962 

0.8543 

0 8229 

0. 7974 

0.9605 

8 

0.9784 

0.9140 

0. 8715 

0.8391 

0.8129 

0.9812 

9 

0.9944 

0.9300 

0.8867 

0 8536 

0.8266 

0.9999 

1G 

i.  0092 

0 . 9445 

0.9005 

0 8666 

0.8390 

1 . 0 1 69 

12 

1.0358 

0.9701 

0 9245 

0. 8893 

0. 8606 

1.0468 

14 

1.0588 

0.9918 

0. 9450 

0.9086 

0 8792 

1.0726 

46 

1-0794 

1.0110 

0.9829 

0. 9255 

n aach 

» • W / ^/ V 

1 AO  C 4 
a.  w 7 

18 

1.0978 

1.0280 

0.9789 

0 9408 

0.9101 

1.  1157 

20 

1. 1143 

1.0436 

0.9931 

0.9546 

0.9235 

i. 1339 

25 

1.  1501 

1.0766 

1.0244 

0.9844 

0.9526 

1.  1737 

30 

1.  1802 

1.  1044 

1.0507 

1.0097 

0.9772 

1. 2070 

35 

1.2060 

1.  1281 

1.0736 

1.0320 

0.9992 

1. 2359 

40 

1 . 2286 

1. 1494 

1.0940 

1 0520 

1 0.191 

1. 2615 

45 

1-2491 

1.  1685 

1 1125 

1.0703 

1.0372 

1 . 2848 

50 

1.2678 

1.  1862 

1. 1296 

1.0871 

1.0538 

1. 3063 

60 

1.3008 

1.2176 

1. 1604 

1. 1179 

1.0844 

1. 3446 

70 

1.3297 

1.2455 

1. 1877 

1 1450 

1.1115 

1 . 3789 

80 

1.3555 

1.7.708 

1.2125 

1 1699 

1. 1363 

1.4095 

90 

1 3789 

1.29  3 8 

1. 2354 

l.  1926 

1. 1595 

1.4377 

100 

1.4009 

1.3150 

1.2568 

1. 2141 

1.  1809 

1.4637 

2®0 

1-5626 

1.4758 

1 4185 

1 3768 

1 3452 

1 6611 
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a = 14 


T* 

ZU-  2) 

Z<2.  3) 

Z<2>  4) 

Z(2,  5) 

ZU.  6) 

’ Z(4’  4) 

0 

1.  1947 

1.0951 

1.0221 

0.9706 

0.9193 

1.0927 

. 1 

1.  2017 

1.  1027 

1.0301 

0.9742 

0.9296 

1.0990 

.2 

1.2090 

1. 1128 

1.0416 

0.9831 

0.9309 

1. 1098 

.3 

1.2109 

1.  1046 

1.0168 

0.9391 

0.8695 

1.0965 

.4 

1. 1947 

1.0716 

0.9690 

0.8825 

0.8109 

1.0607 

.5 

1.  1670 

1.0317 

0.9235 

0.8382 

0.7718 

1.0206 

.6 

1.  1359 

0.9948  ' 

0.8878 

0.8076 

0.7481 

0.9846 

.7 

1.1060  . 

0.9641 

0.8615 

0.7877 

0.7346 

0.9549 

.8 

1.0794 

0.9397 

0 . 8428 

0 . 7754 

0.7277 

0.9315 

.9 

1.0564 

0.9209 

0.8300 

0.7682 

0.  7249 

0.9135 

1.0 

1.0370 

0.9066 

0.8215 

0 . 7647 

0.7250 

0 . 900 1 

1.2 

1.0078 

0.8882 

0.8132 

0.7641 

0.7299 

0.8829 

1.4 

0.9885 

0.8790 

0.8122 

0.7686 

0.7380 

0 . 8746 

1.6 

0.9761 

0.8756 

0.8151 

0 . 7756 

0.7476 

0 . 8720 

1.8 

0.9688 

0 . 8759 

0.8205 

0.7840 

0.7580 

0.8731 

2.0 

0.9649 

0.8785 

0.8271 

0.7930 

0.7683 

0.8764 

2.5 

0.9651 

0.8910 

0.8464 

0.8162 

0.7937 

0.8905 

3.0 

0.9729 

0.9069 

0.8666 

0.3385 

0.8172 

0.9081 

- 3.5 

0.9844 

0 . 9239 

0.8862 

0 . 8594 

0.8384 

0.9267 

> 4 

0.9973 

0.9409 

0.9048 

0.8785 

0 . 8575 

0.9449 

5 

1.0244 

0.9729 

0.9384 

O.f 121 

0.8908 

0 . 9790 

6 

1.0507 

1.0015 

0.9675 

0.9408 

0.9188 

1.0092 

7 

1.0755 

1.0274 

0.9882 

0.9656 

0.9428 

1.0361 

8 

1.0982 

1.0504 

1.0154 

0.9875 

0.9639 

1.0600 

9 

1.  1192 

1.0712 

1.0357 

1.0068 

,0.9827 

1.0816 

10 

1. 1385 

1.0902 

1.0539 

1.0244 

0.9994 

1.  1013 

n 

1. 1731 

1.1237 

1.0859 

1.0554 

1.0293 

1. 1359 

14 

1.2034 

1. 1526 

1. 1136 

1.0818 

1.0550 

1. 1656 

16 

1.2302 

1. 1781 

1.  1379 

1.  1051 

1.0777 

1. 1917 

18 

1.  2540 

l.  2007 

1. 1596 

1. 1261 

1. 0980 

* ^ i r ^ 

l . 

20 

1.2758 

1.2213 

1.  1792 

1. 1452 

1. 1168 

1 . 2364 

25 

1.3229 

1.2659 

1.2222 

1.1868 

1. 1576 

1. 2826 

30 

1.  3622 

1.3034 

1. 2583 

1.2222 

1.  1925 

1.3217 

35 

1.3963 

1.3363 

1.290  2 

1. 2535 

1. 2231 

1. 3559 

1 40 

1.4266 

1.3653 

1. 3185 

1.2816 

1.2510 

1.3863 

45 

1.4543 

1.3917 

1.3445 

1. 3070 

1.2766 

1.4141 

50 

1.4793 

1.4161 

1.3684 

1.  3308 

1 . 2999 

1.4396 

60 

1.5239 

1.4596 

1.4115 

1. 3735 

1. 3427 

1.4856 

70 

1. 5636 

1.4983 

1.4497 

1.4118 

1.3806 

1.5263 

SO 

1 . 5990 

1.  5332 

1.4846 

1.4463 

1.4156 

1,5633 

90 

1.6314 

1. 5656 

1. 5165 

1.4782 

1.4474 

1 . 5972 

100 

1.6612 

1.  5950 

1. 5461 

1.5079 

1.4771 

1.6289 

200 

1.  8858 

1.8198 

i. 7713 

1 7344 

1.7043 

1.8674 
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<X  = 15 


T* 

ZU. 

zd»  2) 

Zd.  3) 

zd.  4> 

Zd.  5) 

z<3*  3> 

0 

1. 1870 

1.0551 

0.9672 

0.9027 

0.8572 

1. 1178 

. 1 

1.  1722 

1.0378 

0,9477 

0.8S10 

0. 8285 

1.  1007 

.2 

1. 1577 

1.0163 

0.9153 

0.8343 

0.7655 

1.0765 

.3 

1.  1320 

0.9750 

0.8594 

0.7687 

0.6970 

1.0392 

.4 

1.0984 

0.9300 

0.8099 

0.7221 

0.6579 

1.0012 

.5 

1.0645 

0.8916 

0.7747 

0 . 6944 

0.6388 

0.9672 

.6 

1.0337 

0.8618 

0.7514 

0. 6792 

0.6396 

0.9387 

.7 

1.0075 

0.8397 

0.7367 

0.6715 

0.6285 

0.9158 

.8 

0.9857 

0.8235 

0.7278 

0.6687 

0.6299 

0.8978 

9 

0.9678 

0.8121 

0.7230 

0.6687 

0.6333 

0.8840 

1.0 

0.9533 

0 . 8040 

0. 7208 

0.6707 

0.6378 

0.8735 

1.2 

0.9323 

0.7954 

0.7218 

0.6777 

0 . 6484 

0.8605 

1.4 

0.9191 

0.7931 

0.7266 

0.6867 

0.6596 

0 . 8548 

1.6 

0.9112 

0 . 7944 

0.7335 

0.6964 

0.6708 

0.8535 

1.8 

0.9070 

0.7979 

0.7412 

0.7063 

0.6817 

0.8554 

2.0 

0.9054. 

0.8028 

0 . 7494 

0.7159 

0.6921 

0. 8593 

2.  5 

0.9083 

0.8179 

0.7699 

0.7389 

0.7160 

0 . 8734 

3.0 

0.9167 

0.8344 

0. 7896 

0.7597 

0.7372 

0.8903 

3.5 

0.9273 

0.8507 

0.8078 

0.7783 

0.7558 

0.9077 

4 

0.9389 

0.8662 

0.8245 

0.7953 

0.7725 

0.9245 

5 

0.9627 

0.8951 

0.8541 

0.8245 

0. 8012 

0.9558 

6 

0.9856 

0.9203 

0.8796 

0 . 8494 

0.8250 

0.9836 

7 

1.0069 

0.9430 

0.9017 

0.8708 

0.8458 

1.0086 

8 

1.0268 

0.9632 

0.9214 

0.8895 

0.8638 

1.0308 

9 

1.0449 

0.9816 

0.9390 

0.9063 

0.8800 

1.0514 

10 

1.0619 

0.9981 

0.9549 

0.9216 

0.8946 

1.0699 

12 

1.0922 

1.0277 

0.9829 

0.9483 

0.9205 

1.  1032 

14 

1. 1188 

1.0530 

1.0070 

0.9715 

0.9428 

1. 1321 

16 

1. 1426 

1.0756 

1.0284 

0.9921 

0.9628 

1. 1579 

18 

• • / * *. 

1.  1QW 

1. 0957 

1 . 0475 

1.0104 

0.9805 

1. 1812 

20 

1. 1834 

1. 1140 

1.0648 

1.0273 

0.9970 

1.  20  26 

25 

1.2257 

1. 1537 

1.1029 

1.0642 

1.0333 

1. 2494 

30 

1.2614 

1. 1872 

1.  1353 

1.0957 

1.0642 

1. 2897 

35 

1.  29  23 

1. 2167 

1.  1637 

1. 1234 

1.0920 

1 . 3254 

40 

1. 3199 

1 . 2429 

1.1890 

1. 1486 

1 .1169 

1. 3574 

45 

1. 3446 

1. 2666 

1. 2123 

1. 1718 

1.  1397 

1. 3869 

50 

1. 3675 

1. 2886 

1. 2339 

1.  1929 

1.  1609 

1.4141 

60 

1.4083 

1 . 3280 

1.2726 

1.2316 

1.  1996 

1.4631 

70 

1.4443 

1. 3630 

1. 3073 

1. 2660 

1. 2340 

1.5072 

80 

1.4766 

1. 39  45 

i. 3389 

1.2974 

1. 2654 

1.5471 

90 

1.50  59 

1.4235 

1. 3675 

1.3267 

1 . 294b 

1. 5837 

100 

1.5331 

1.4505 

1.  3945 

1.3534 

i.  3217 

1. 6180 

200 

1. 7372 

i 

1. 6533 

1. 5991 

1. 5596 

1 . 529 1 

1.  8773 

Table  IV. 


Collision  Integral  Functions  continued 


WIS-ONR-  1 
19  June  1953 


38 


OL  = 15 


jr' 

t 


T* 

z<2,  2) 

z(2,  3) 

! ZU-  4) 

z<2.  5) 

z(2,  6) 

ZU,  4) 

0 

1.  1947 

1 0951 

1.0221 

0.9706 

0.9193 

1.0927 

. 1 

1 2016 

1.  1025 

1.0306 

0.9736 

0.9290 

1 1061 

. 2 

1.2081 

1 1114 

1.0403 

0.9828 

0 9325 

1.  1150 

. 3 

1.  2107 

1.  1061 

1-0211 

0.9465 

0.8796 

i in  ->n 

+■  • M / 

. 4 

1.  1978 

1.0783 

0.9791 

0.8952 

0.8252 

1.0705 

. 5 

1.  1739 

1 0426 

0 9374 

0.8536 

0.7879 

1 0333 

. 6 

1. 1462 

1.0088 

0.9039 

0.8247 

0. 7654 

0 9996 

. 7 

1. 1191 

0.9803 

0 879  1 

0 8060 

0 7529 

0.9718 

. 8 

1.0945 

0 9575 

0.8617 

0. 7946 

0 7470 

0.9499 

.9 

1.0734 

0.94C0 

0.8498 

0.7883 

0 . 7454 

0.9331 

1.0 

1.0555 

0 9268 

0.8423 

0.7858 

0.7466 

0.9209 

1.2 

1.0285 

0 9102 

0.8359 

0. 7872 

0.7535 

0.9057 

1.4 

1.0110 

0 9027 

0.8365 

0.7936 

0. 7637 

0.8992 

1.  6 

1.0003 

0.9009 

0.8412 

0.8024 

0 . 7751 

0. 8982 

1.  8 

0.9944 

0 9027 

0,8481 

0 8124 

0.7868 

0.9007 

2.0 

0.9920 

0.9067 

0 8562 

0 8227 

0.7986 

0 9054 

2.  5 

0 9953 

0 9223 

0.8786 

0.8488 

0.8266 

0.9222 

3.0 

1.0059 

0 9409 

0.9013 

0 8734 

0.8521 

0.9419 

3.  5 

1.0196 

0.960  2 

0.9229 

0.8960 

0 8750 

0.9622 

4 

1.0345 

0 9790 

0.9431 

0 9169 

0. 8958 

0.9819 

5 

1.0651 

1 0139 

0,9796 

0 9533 

0.9318 

1.0187 

6 

1.0943 

1,0452 

1.0111 

0.9844 

0.9623 

1.0514 

7 

1. 1213 

1,0733 

1 0389 

1.0115 

0 9887 

1 0805 

8 

1. 1462 

1.0985 

1.0634 

1.0356 

1.0121 

1.  1064 

9 

1. 1691 

1. 1212 

1.0856 

1.0570 

1.0332 

1. 1300 

10 

1. 1904 

1 1421 

1.  10  59 

1 0766 

1.0523 

1. 1515 

12 

1. 2283 

1. 1791 

1. 1416 

1. 1115 

1.0860 

1.  1895 

14 

1.2616 

1 2111 

1. 1727 

1. 1416 

1,1154 

I 2225 

i a 

A.  w 

1 90  11 

* • • 

1.  239? 

1 200  2 

1. 1683 

1 » 14) 7 

1 . 2518 

18 

1. 3179 

1. 2655 

1.2250 

1.1925 

1. 1656 

1 2781 

20 

1 . 3423 

1. 2887 

1.2477 

1.2148 

1 1873 

1. 3020 

25 

1.3954 

1 3399 

1.2975 

1 2634 

1. 2353 

1. 3544 

30 

1.4406 

1 3835 

1.3400 

1. 30  54 

1. 2768 

1.  3992 

35 

1 . 480  2 

1. 4219 

1 3775 

1, 3425 

1. 3138 

1.4382 

40 

1. 5152 

1.4561 

1 4114 

1 3759 

1. 3471 

1.4737 

45 

1. 5474 

1. 4874 

1 4424 

1.4067 

1. 3776 

1. 5057 

50 

1.5769 

1 5163 

1.4710 

1 4352 

1.4060 

1. 5353 

60 

1. 6298 

I. 5684 

1. 5225 

1. 4865 

1.4575 

1. 5889 

70 

1. 6768 

1.6150 

1,5687 

1 5326 

1. 5034 

1 6365 

80 

I.  7190 

1.6569 

1.6107 

1 5747 

1.5453 

j / 

i , 017  1 

90 

l.  7580 

1 6952 

1.6491 

1.6132 

1.5841 

1.  7194 

100 

1. 7938 

1 7312 

1. 6850 

1, 6489 

1.6200 

1.  7562 

200 

2 0647 

2 0026 

1 9572 

1.9227 

1.8945 

2.0351 

WIS-ONR-i 

Table  V Functions  for  calculating  the  higher  approximations  to  19  June  1953 

the  transport  coefficients  . 39 


The 

function  f 

T* 

cl  = 12 

OL  = 13 

«=  14 

°L  = 15 

0 

l 0017 

1 0017 

1 0017 

1.0017 

. 1 

1.0019 

1 0018 

1.0018 

1 0018 

. 2 

1 00  17 

1.00  21 

1.0021 

1.0020 

. 3 

1 0012 

1.00 14 

1 0015 

1 0016 

.4 

1 000  3 

1 0005 

1.0006 

1.0008 

5 

1.0001 

1.0001 

1 .0002 

1.0003 

.6 

1.0001 

1.0000 

1.0000 

1.0001 

. 7 

1.000  2 

1.0001 

1 0000 

1.0000 

. 8 

1.0002 

1 0001 

1 0000 

1.0000 

.9 

1.0002 

1.000  1 

1.0000 

1.0000 

1.0 

1 . 000  1 

1.0001 

1. 0000 

1 0000 

1.2 

1 0001 

1 000  1 

1.0C01 

1.0001 

1.4 

1 0001 

1.0002 

1 0002 

1.0003 

1.6 

1.0003 

1.0004 

1.0005 

1.0007 

1.8 

1 000  5 

1 0007 

1.0008 

1.0010 

2,0 

1.0008 

1.001 1 

1.0012 

1.0015 

2.  5 

1.0017 

1 0020 

1.0022 

1.0026 

3.0 

1 00  25 

1.0029 

1.0031 

1.0035 

■ 3.5 

1.0032 

1.0036 

1.0039 

1,0043 

4 

1.0037 

1.0042 

1.0046 

1.0050 

5 

1.0045 

1 0050 

1. 0055 

1.0058 

6 

1 0050 

1.0055 

1.0061 

1.0064 

7 

1.0052 

1 nn  e 0 

1 . V/W-/U 

1 0066 

1 0067 

8 

1 00  54 

1.0060 

1 0067 

1 0070 

9 

1.0055 

1.0061 

1.0068 

1.0071 

10 

1 0056 

1 0062 

1 0069 

1.0072 

12 

1 00  57 

1.0062 

1.0069 

1 0073 

14 

1 0057 

1.0063 

1.0069 

1 0073 

16 

1.0057 

1 0063 

1 0069 

1.0073 

18 

1 nn  K7 

* • v w 1 

1 nn  L a 

X W 

i nn  L n 
X c WVU7 

1 AAHO 

A . UU  | J 

20 

1.0057 

1.0064 

1.0068 

1.0073 

25 

1 0057 

1.0064 

1. 0068 

1.0073 

3€ 

1.0057 

1. 0065 

1.0068 

1. 0074 

35 

1.0057 

1.0065 

1. 0068 

1.0074 

40 

1.0058 

1.0066 

1 0068 

1.0075 

45 

1 0058 

1 0067 

1.0068 

1.0075 

50 

1.0059 

1.0068 

1.0068 

1.0075 

60 

1.0060 

1.0069 

1.0069 

1.0077 

70 

1.0061 

1 0071 

1.0070 

1.0078 

80 

1 0062 

1.0072 

1.0071 

i . 0079 

90 

1 0063 

1 0073 

1 0072 

1. 0080 

100 

1.00  64 

1 0074 

1. 0073 

1. 0082 

200 

1 0074 

1. 0084 

1.0081 

1. 0091 

Table  V. 


Functions  for  calculating  the  higher  approximations  to 


WIS-ONR-  1 
19  June  1953 


the  transport  coefficients,  continued 


The 

function  f^ 

T* 

= 12 

= 13 

= 14 

= 15 

0 

1.0027 

1. 0027 

1.0027 

1 0027 

. 1 

1 0029 

1. 0028 

1.0029 

1.0028 

. 2 

1 0027 

1.00  32 

1 0033 

1.0032 

. 3 

1.0020 

1.0023 

1.0025 

1.0026 

.4 

1.0006 

1.0009 

1,0011 

1.0013 

. 5 

1.0002 

1.0002 

1.0003 

1.0005 

. 6 

1.0002 

1 . 000  1 

1.0001 

1.0001 

. 7 

1.0003 

1.0001 

1.0000 

1.0000 

. 8 

1. 0003 

1.0001 

1.0000 

1.0000 

.9 

1.0003 

1 0001 

1.0000 

1.0000 

1.0 

1.0002 

1.0001 

1.0000 

1.0000 

1.2 

1.0002 

1 . 000  1 

1.0001 

1.0002 

1.4 

1. 0002 

1.0003 

1 0004 

1.0005 

1.6 

1.0004 

1.0006 

1.0008 

1.0010 

1.8 

1.0008 

1.0011 

1.0013 

1.0016 

2.0 

1.0013 

1.0016 

1.0019 

1.0022 

2.  5 

1.0026 

1 0031 

1.0034 

1.0039 

3.0 

1.0038 

1.0045 

1.0048 

1.0054 

3.  5 

1.0049 

1. 0056 

1.0060 

1.0067 

4 

1.0058 

1.0065 

1.0071 

1.0077 

5 

1.0070 

1.0078 

1.0085 

1.0091 

6 

1 0077 

1. 0085 

1.0094 

1.0099 

7 

1.0081 

1 0090 

1 0105 

1.010  5 

8 

1.0085 

1.0093 

1.0104 

1.0109 

9 

1.0087 

1.0095 

1.0106 

1.0111 

10 

1.0088 

1 0097 

1.0107 

1.0112 

12 

1.0089 

1.0098 

1. 0109 

1.0114 

14 

1.0089 

1 0099 

1 0108 

1.0114 

16 

1.0089 

1.0099 

1.0108 

1.0115 

1 ft 
A O 

1 . GOSS 

1. 0100 

1.0108 

1.0115 

20 

1.0089 

1.0100 

1.0107 

1.0114 

25 

1.0089 

1.0101 

1.0106 

1.0115 

30 

1.0089 

1.0101 

1.0106 

1.0115 

35 

1.0090 

1.0102 

1. 0106 

1.0116 

40 

1 0090 

1.0103 

i.0106 

1.0116 

45 

1.0091 

1.0104 

1. 0106 

1.0117 

50 

1.0092 

1. 0105 

1.0107 

1. 0118 

60 

1.0093 

1.0108 

1.0108 

1. 0120 

70 

1.0096 

1.0111 

1. 0109 

1.0122 

80 

1.0097 

1 0112 

1.0110 

1.0124 

90 

1.0099 

1 0114 

1.0113 

1.0125 

100 

1.0101 

1.0116 

1.0114 

1 0127 

200 

1 0115 

1.0132 

1 0127 

1 0141 

WIS-ONR-  1 
19  June  1953 

Table  V.  Functions  for  calculating  the  higher  approximations  to 


the  transport  coefficients,  continued 


The  function 

T* 

= 12 

= 13 

= 14 

= 15 

0 

1.0020 

1.0020 

1.0020 

1.0020 

. 1 

1.0017 

1.0017 

1.0018 

1.0018 

.2 

1.0010 

1.0012 

1.0013 

1,0014 

. 3 

1.0002 

1.0004 

1.0005 

1.0006 

.4 

1.0001 

1.0001 

1.0001 

1.0002 

. 5 

1.0002 

1.0001 

1.0000 

1.0000 

.6 

1.0003 

1.0001 

1.0000 

1.0000 

.7 

1.0003 

1.0001 

1.0000 

1.0000 

.8 

1.0002 

1.0001 

1.0000 

1.0000 

.9 

1.0002 

1.0001 

1.0000 

1.0000 

1.0 

1.0001 

1.0000 

1,0000 

1.0001 

1.2 

1.0001 

1.000  1 

1.0002 

1.0003 

1.4 

1.0001 

1.0003 

1.0004 

1.0006 

1.6 

1.0003 

1.0005 

1.0007 

1.0009 

1.8 

1.0006 

1.0008 

1.0011 

1.0013 

2.0 

1.0009 

1.0012 

1.0015 

1.0017 

2.5 

1.0017 

1.0020 

1 . 00  24 

1.0028 

3.0 

1.0024 

1.0028 

1.0032 

1.0037 

3.  5 

1.0030 

1.0035 

1.0039 

1.0044 

4 

1.0035 

1.0041 

1.0045 

1.0050 

5 

1.0042 

1.0049 

1.0053 

1.0059 

6 

1.0047 

1.0053 

1.0059 

1.0064 

7 

1.0049 

1.0057 

1.0062 

1.0068 

8 

1.0051 

1.0059 

1.0065 

1.0070 

9 

1.0052 

1.0060 

1.0066 

1.0072 

10 

1.0053 

1.0061 

1.0067 

1.0073 

12 

1.0053 

1.0061 

1.0068 

1.0074 

14 

1.0053 

1 0062 

1.0069 

1.C075 

16 

1.0052 

1.0061 

1.0068 

1.0075 

IS 

1.0052 

1.0061 

1.0068 

1.0075 

20 

1.0051 

1.0060 

1.0068 

1.0075 

25 

1.0050 

1.0059 

1.0068 

1.0075 

30 

1.0049 

1.0058 

1.0067 

1.0075 

■'  35 

1.0048 

1.0058 

1.0066 

1.0075 

40 

1.0048 

1.0058 

1.0067 

1.0075 

45 

1.0047 

1.0058 

1.0066 

1.0075 

50 

1.0047 

1.0058 

1. 0067 

1.0076 

60 

1.0047 

1.0058 

1.0067 

1.0077 

70 

1.0047 

1.0058 

1.0068 

1.0077 

80 

1.0047 

1.0059 

1.0069 

1.0078 

90 

1.0047 

1.0059 

1.0070 

1.0079 

100 

1.0048 

1.0060 

1.0070 

1.0081 

200 

1.0053 

1.0067 

1.0077 

1 

1.  0088 

Table  VI. 


The  reduced  thermal  diffusion  ratio 


W1S-ONR-  1 
19  June  1953 


9 


for  calculating  isotopic  thermal  diffusion. 


T T 


M,  - Ma. 

— — ^ X x 

M,  +M,t 


42 


T* 

oc-12 

c:«13 

oc-14 

a = 15 

0 

0. 3133 

0.3133 

0. 3133 

0 3133 

. 1 

0. 2844 

0 2865 

0 . 2887 

0. 2887 

. 2 

0. 2047 

0 2264 

0 . 240  1 

0. 2468 

. 3 

0.0738 

0, 1110 

0.  1369 

0. 1545 

.4 

-0.0271 

0.0157 

0.0484 

0.0727 

5 

-0.0829 

-0  0389 

-0.0037 

0.0229 

.6 

-0.  1038 

-0.0616 

-0.0256 

0.0020 

. 7 

-0. 1031 

-0.0623 

-0.0270 

0.0005 

.8 

-0.0902 

-0 .0  509 

-0.0162 

0.0105 

-9 

-0.0690 

-0  0317 

0.0023 

0.0  287 

1.0 

-0.0456 

-0.0088 

0.0244 

0.0498 

1.  2 

0.0071 

0.0419 

0.0734 

0 0981 

1.4 

0.0580 

0 0923 

0. 1226 

0.  1466 

1.  6 

0.  1050 

0.1386 

0.  1680 

0.  1911 

1.8 

0 1476 

0. 1797 

0. 2083 

0 2310 

2.0 

0. 1841 

0.2161 

0.  2438 

0.2665 

2.  5 

0. 2570 

0.2883 

0.3145 

0 3380 

3.0 

0. 3103 

0 3404 

0.3658 

0.3890 

3.  5 

0 . 3469 

0 3786 

0. 4035 

0.4270 

4 

0.3744 

0.4060 

0.4313 

0 4544 

5 

0.4111 

0.4437 

0.4679 

0.4931 

6 

0.4312 

0. 4643 

0.4911 

0. 5134 

7 

0.4414 

0.4778 

0. 5044 

0. 5283 

8 

0 . 4494 

0.4358 

0.5131 

0. 5358 

9 

0.4531 

0.4889 

0.  5181 

0.  5427 

10 

0.4536 

0.4911 

0. 5208 

0.  5446 

12 

0.4536 

0 4917 

0.  5233 

0. 5494 

14 

0. 4512 

0.4919 

0. 5225 

0 5495 

16 

0 4469 

0.4887 

0. 5209 

0. 5496 

18 

0 . 4441 

0.4866 

0. 5187 

0.  5485 

20 

0.4387 

0.4826 

0 . 5190 

0 5479 

25 

0. 432.1 

0 4764 

0 . 5i40 

0. 5452 

30 

0 . 4245 

0 4713 

0. 5106 

0. 5427 

35 

0 4188 

0.4681 

0. 5066 

0. 5429 

40 

0.4161 

0 4648 

0 50  58 

0 . 5426 

45 

0.4118 

0 4641 

0.  50  38 

0.  5426 

50 

0 409  3 

0.4616 

0. 50  35 

0 5429 

60 

0 4051 

0.4598 

0. 5032 

0. 5443 

70 

0. 4048 

0 4590 

0. 5044 

0 . 5460 

80 

0.4034 

0.4583 

0 5071 

0 . 5474 

90 

0 4022 

0.4584 

0. 5093 

0.  550  5 

100 

0 . 4024 

0.4586 

0 5095 

0. 5534 

200 

0 4110 

0 . 4703 

0.  5277 

0.  5727 

Table  Vll  Functions  for  calculating  the  transport  properties 


WIS-ONR-  1 
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of  mature  s 


43 


The  function  A 

T+ 

— 

u = 12 

ot  _ 13 

ot  = 14 

<*=  15 

0 

1.0065 

1.0065 

1 0065 

1.0065 

. 1 

1 0060 

1 0146 

1 0 247 

1.0251 

*> 

<*>  1 

1.0  374 

1 0374 

1 . 0445 

1.0435 

3 

1 0731 

1 .0693 

1 0729 

1 0695 

.4 

1.0976 

1,0932 

1.0946 

1 090  5 

.5 

1 1108 

1.  1067 

1. 1070 

1 1028 

6 

1 1170 

1.1131 

1.1128 

1.  1088 

7 

1 1193 

1.1154 

1 1145 

1.  1108 

.8 

1 1194 

1.1154 

1.  1142 

1 1104 

.9 

1. 1185 

1 1143 

1 1126 

1.  1091 

1.0 

1. 1173 

1. 1126 

1 1106 

1.  1071 

1 2 

1.  1146 

1 1096 

1 1067 

1.  1032 

1.4 

1.  1126 

1 1070 

1.  1037 

1.  1000 

1.6 

1.1114 

1.  1052 

1.1013 

1.0978 

1 8 

1.1111 

1. 1044 

1 0999 

1.0964 

2.0 

1 1111 

1.  1044 

1.0990 

1 0956 

2 5 

1 1130 

1 1056 

1, 0987 

1 0958 

3 0 

1 1161 

1 1081 

1.  1002 

1.0973 

3 5 

1.1193 

1.1113 

1. 1026 

1.0995 

4 

1. 1226 

1 1141 

1 1050 

1 1018 

5 

1. 1284 

1 1 194 

1 1101 

1. 1064 

6 

1.  1335 

1 1238 

1 1147 

1.  1103 

7 

1. 1377 

1. 1277 

1.  1189 

1 . 1 136 

8 

1.  1417 

1 1309 

1 1224 

1 1163 

9 

1 1449 

1 1336 

1. 1255 

1 1189 

10 

1 1477 

1 1361 

1 1281 

1. 1210 

12 

1 1532 

1 1406 

1 1326 

1 1246 

14 

1.  15&3 

1 1447 

1 1366 

1 1276 

16 

1.1625 

1 1485 

1 1397 

1 1300 

18 

l 1665 

1 . 1517 

1 1423 

1. 1322 

20 

1.  1702 

1 1548 

1 1449 

1. 1343 

25 

1 1790 

1 1623 

1 150  2 

1 1385 

30 

1.-1865 

1. 169  1 

1.  1542 

1.  1421 

35 

1 1938 

1 1752 

1 1578 

1 1454 

40 

1 200  5 

i 3813 

1. 1612 

1. 1480 

45 

1 2067 

1 1867 

1 1643 

1 . 1508 

50 

1 2123 

1 1917 

1 1668 

1. 1531 

60 

1 2230 

1 200  8 

1 1715 

1 1573 

70 

1 2^26 

1 209  2 

1 1759 

1 1610 

80 

1 2416 

1 2168 

1 1796 

1 1642 

90 

1 2492 

1 2240 

1.  1831 

1 1674 

100 

1 2569 

1 2300 

1 1858 

1 1700 

200 

1.  3103 

1 2754 

i . 2068 

1.  1885 

Table  VII.  Functions  lor  calculating  the  transport  properties 


WLS-ONR-  1 
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of  mixtures  , continued. 


44 


1 

Tne  junction  a 

1 

T* 

<*  = 12 

ot  = 13 

d=  14 

<*=  15 

0 

1. 1851 

1.  1851 

1.  1851 

1. 1851 

. 1 

1.  1988 

1.  1955 

1.  1926 

1.  19  28 

. 2 

1 . 2528 

1. 2407 

1. 2328 

1.  2268 

. 3 

1. 3028 

1. 2895 

1. 2788 

l.  2698 

.4 

1.  3127 

1. 3019 

1.2922 

1.  2840 

. 5 

1. 300  3 

1. 2922 

1.  2844 

1. 2768 

.6 

1. 2808 

1.2737 

1. 2668 

1. 2609 

.7 

1. 2601 

1.2539 

1.2476 

1 . 2424 

. 8 

1.2403 

1. 2346 

1. 2289 

1. 2238 

.9 

1.2239 

1. 2172 

1. 2125 

1. 2074 

1.0 

1. 2079 

1. 20  2 5 

1.  1973 

1.  1930 

1.  2 

1. 1837 

1. 1780 

1. 1733 

1.  1689 

. 1.4 

1. 1662 

1. 1612 

1. 1561 

1. 1523 

1.6 

1.1536 

1. 1490 

1.  1440 

1.  1392 

1.8 

1. 1458 

1.  1394 

1.  1347 

1.  1298 

2.0 

1. 1387 

1.  1324 

1.  1279 

1. 1226 

2.5 

1.1288 

1.  1222 

1.  1174 

1. 1119 

3.0 

1-  1250 

1.  1167 

1.1113 

1.  1057 

3.  5 

1. 1217 

1.  1136 

1.  1084 

1.  10  24 

4 

1. 1211 

1.  1128 

1.  1067 

1.  1002 

5 

1. 1225 

1.  1126 

1.  1048 

1.  1001 

6 

1. 1242 

1.  1131 

1.  1059 

1 . 0939 

7 

1. 1252 

1.  Ii55 

1.  1067 

l.  1006 

8 

1. 1278 

1. 1176 

1. 1079 

1.  1009 

9 

1. 1298 

1.  1185 

1.  1094 

1.  1025 

10 

1. 1317 

1 1203 

1.  1103 

1.  1026 

12 

1. 1351 

1.  1225 

1. 1127 

1.  1050 

14 

1 1376 

1. 1247 

1.  1135 

1.  1056 

16 

1. 1384 

1. 1259 

1.  1149 

1. 1066 

18 

1.  1405 

1. 1271 

1. 1153 

1.  1070 

20 

1. 1404 

1. 1274 

1. 1178 

1.  1077 

25 

1.  1426 

1. 1282 

1. 1176 

1. 1070 

30 

1. 1426 

1. 1281 

1. 1178 

1.  1058 

35 

1 . 1427 

1. 1279 

1. 1162 

1.  1055 

40 

1. 1427 

1.  1268 

1.1159 

1. 1050 

45 

1. 1418 

1.  1270 

1. 1148 

1.  1035 

50 

1.  1407 

1. 1255 

1.  1142 

1.  10  23 

60 

1 . 1388 

1.  1236 

1.1119 

1.  1003 

70 

1. 1378 

1. 1218 

l.  1106 

1.0980 

80 

1. 1365 

1.1198 

1.  1096 

1.0950 

90 

1.1350 

1. 1180 

1.  1077 

1.0940 

100 

1. 1331 

1. 1154 

1.  1049 

1.0922 

200 

1.  1206 

1.  1021 

1.  0911 

1. 0 779 

Table  VU*  Functions  for  calculating  the  transport  properties 
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The  function  C* 

T* 

M 

II 

a.  = 13 

<*=  14 

a=  15 

0 

0.8889 

0.8889 

0.8889 

0.8889 

. 1 

0.8335 

0.8844 

0.8853 

0.8853 

.2 

0.8697 

0.8736 

0.8765 

0.8779 

.3 

0.8463 

0.8531 

0.8581 

0.8613 

.4 

0.8279 

0.8360 

0.8421 

0.8467 

.5 

0.8172 

0.8257 

0.8325 

0.8376 

.6 

0.8128 

0.8212 

0.8283 

0.833? 

.7 

0.8127 

0.8209 

0 . 8280 

0.8334 

.8 

0.8151 

0.8231 

0.8301 

0.8354 

.9 

0.8193 

0.8269 

0.8338 

0.8391 

j 1.0 

0.8240 

0.8315 

0.8383 

0.8434 

•\  1.2 

0.8348 

0.8419 

0.8482 

0,8532 

1.4 

0.8453 

0.8521 

0.8582 

0.8629 

1.6 

0.8549 

0.8615 

0.8673 

0.8718 

1.8 

0.8636 

0.8698 

0.8753 

0.8797 

2.0 

0.8709 

0.8771 

0.8823 

0.8867 

2.5 

0.8855 

0.8913 

0.8961 

0.9005 

3.0 

0.8959 

0.9014 

0.9059 

0.9102 

3.5 

0.9031 

0.9088 

0.9131 

0.9174 

4 

0.9084 

0.9141 

0.9183 

0.9226 

5 

0.9155 

0.9213 

0.9253 

0.9298 

6 

0.9195 

0.9254 

0.9298 

0.9337 

7 

0.9216 

0.9280 

0.9324 

0.9365 

8 

0.9232 

0.9295 

0.9342 

0.9381 

9 

0.9241 

0.9302 

0.9352 

0.9394 

10 

0.9242 

0.9307 

0.9359 

0.9399 

12 

0.9245 

0.9311 

0.9366 

0.9409 

14 

0.9243 

0.9313 

0,9367 

0.9412 

L6 

ft  09**7 

n o*i a 

w • / y * v 

A 012.2 
V • 7 rfWW 

A A A \ A 

W « 7-*  ATI 

18 

0.9234 

0.9307 

0.9364 

0.9413 

20 

0.9226 

0.9302 

0.9366 

0.9414 

25 

0.9218 

0.9295 

0.9361 

0.9413 

30  . 

0.9209 

0.9291 

0.9358 

0.9412 

35 

0.9203 

0.9289 

0.9354 

0.9415 

40 

0.9202 

0.9288 

0.9355 

0.9417 

45 

0.9198 

0.9290 

0.9355 

0.9420 

50 

0.9197 

0.9290 

0.9356 

0.9423 

60 

0.9196 

0,9293 

0.9360 

0.9430 

70 

0.9202 

0.9298 

0.9367 

0.9437 

to 

0.9205 

0.9303 

0.9375 

0.9444 

90 

0.9208 

0.9309 

0.9383 

0.9453 

196 

0.9214 

0.9316 

0.9387 

0.9461 

200 

0.9271 

0.9377 

0. 9445 

0. 9520 

Index  for  n and  A 

L 

12  13  ! 14  15 

0 0 833  0 833  0 833  0.833  | 

.1  0.827  0.831  0.830  0.830  j 

2 0 832  0 819  0 818  0.820  1 

. 3 0 874  0 858  if;  0,  851  0 846 

. 4 0 950  0 925  , b . 9 1 2 0 399 

5 1 007  0.980.,  0.964  0 947 

.6  1.040  1.01.4  0.997  0.979 

. 7 1.056  1 03£  1 013  Q.996 

.8  1.058  ;1.03]6'  /:  1.018  1.001 

,9  1.052  1.031  1.013  0.997 

1.0  1.040  1 0 20 ^ 1.003  0 988 

1.2  1.009  p.99l  I 0 975  0 960 

1.4  0 974  p 95%  ‘ 0 943  0 928 

1.6  0.941  p.9Z 0.912  0 897 

1.8  0.911  0.894.''  0,884  0.869 

2.0  0.884  b.86jb;;'  0 858  0 844 

2.5  0.832  P 816",  0^807  0.793 

3.0  0.796  p 78'p 0.771  0.758 

3.5  0.770  j).  75$.'1.  0.746  0.733 

4 0.753  0 73§' J * Oi 726  0.715 

5 0.730  p 715  0.701  0.692 

6 0. 717  0 703  0 687  0. 679 

7 0. 710  0.695'  ; 0.678  0 671 

8 0 706  69<}  0.674  0 666 

9 0 70.3  0 68b(  ‘ 0 672  0 664 

10  0.701  686*  0 670  0.662 

12  0 700  0. 68k'  ‘ 0.669  0.660 

Q (1 i 1 ■ (t 

14  0 700  ) . 683  , 0,.  669  0.660 

16  0.700  P-683!,  0.669  0 659 

18  0.700  682’  0.670  0.659 

20  0 700  0,68%  0.671  0.660 

25  0.700  i) . 680  0.672  0 659 

30  0.699  67$'J  ,0.673  0.659 

35  0.698  67€|  ’ 0.672  0 658 

40  0 697  ) 67& 1 0.672  0.656 

45  0 696  0 674  0.672  0.655 

50  0 694  p 672  V-  0.671  0.654 

60  0.692  t 66$  \ j 0.669  0.651 

70  0.687  P.664  0.667  0.647 

; 80  0 685  0 661  0.665  0.645 

90  0 682  0 659  0,661  0 643 

100  0 679  0 655  0.659  0.640 

200  0 657  0 633  0.640  j 0.620 


Index  for  .1%  ^ 

12  13  M .14  , 'i5 

0 . 833  0 8 33I  ' " P . 8 3 3% ~ 0 , 8 3 3 : 

0.850  ; 0 847;  g?§|40  0!.  844 

0 891  0 8791  0f  87b,.  0.  866 

0 961%,  b..94i:  0-926,.  oj.91% 

1.016  0..992'  0.9*4,  0;.960.  ' 

1 049  1 0 231  ivQQ3()  0 . 987r 

1 062  1,036  0,999, 

1.062  i,-037  ^000\ 

1.0  55  1,  031!  i :,9 10  C)  0,99% 

1.042  1.019i  QpiPo  0 . 9 8 3p 

1.028  1,005  Pt9§5(V  0, 970f) 

0.996  0.974  0,t955$)  0.94%  j 

0 964  0.944!  ?.,^50  0 91%.  J 

0.935  0.915  §t$pn  OL  885^  j 

0 909  *1.891  0T874;,  0,86%  \ 

0.887  0.869  01,  840r>  , 

0.844  0 826  0*^2  ()  0L 799e ’ " 

0 812  0.796  0.  769,/, 

0.791  0,  774|  (j  Ijfei  0 0 . 748^  , 

0.775  0 758  Q ':i  0f73Z  L 

0.753  0.736  0v^4(i  0,711  f 

0.742  0.724!  4^1  o 699 

0 735  0 716|  0,690  ? 

0 730  0.71l!  0(6^7(;  Ot  686  H 

0.728  0,  709;  qt|94  n 0'f  682  (> 

0.727  0.708'  0^69^0  °‘  680  01 

0.727  0,  707  0-677  u 

0.727  0.  7061  4^690  0 0,676  M 

0 729  0,  707  n 676  %, 

0.730  0.708j  (j.^1  '0  0).  676 

0.732  0.710  0.694  0 0,676  ^ 

0 735  0,711  6,j6^2  0 0 676 

0.737  0^713  Q,  fyi  q o[  676  0T 

0.739  713.  o 0,676  u 

0 740  0,714;  o °l675  Oi- 

0.741  0 j 713  QfC&*4  o 0f  674 

0 .741  0!.  713!  0,v69>3  0 0 673  (U. 

0 741  0 , 7 1 2|  O^tftZ  o 0 671  ih> 

0 739  0 [ 710*  0.69,0  -!  01  669  ,.•! 

0 739  0.709  0.  688  , 0.6b7  (;S 

0 738  0.707  0,685  , 0,664  Jv 

0 736  0.705  0 68^  0.662  J0I 

0 719  0.687  0.667  Oi  644  OS 


i 

i 

t 

I 
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Table  VTTTR.  Sutherland  constants  for  viscosity,  thermal  conductivity, 


and  self -diffusion 


S/'i  for  D. 


0.236 
0.231 
0. 229 
0.227 
0. 224 
0.224 
0.222 
0.221 
0. 220 
0. 218 
0.216 
0. 214 
0. 211 
0. 20  8 
0.  202 
0.  197 
0.  192 
0.  189 
0.184 
0.  154 


0. 500 
0.491 
0.467 
0.541 
0. 701 
0.865 
0.988 
1.054 
1.073 
1.054 
1.012 
0.904 
0.796 
0.700 
0.622 
C.558 
0.443 
0.372 
0.326 
0.292 
0.252 
0.230 
0 . 218 
0. 211 
0. 207 
0. 204 
0. 20  3 
0. 20  3 
0. 204 
0. 20  5 
0. 20  6 
0. 208 
0.209 
0. 20  8 
0. 208 
0. 20  8 
0.206 
0.203 
0.  200 
0.  197 
0.  192 
0.  190 
0.  163 


0. 500 
0.492 
0.471 
0.  528 
0 . 664 
0.810 
0.921 
0.985 
1.003 
0.989 
0.952 
0.852 
0.750 
0 . 660 
0. 584 
0.524 
0.415 
0.349 
0. 304 
0.273 
0. 238 
0.219 
0. 207 
0. 200 
0.  196 
0.  194 
0.  191 
0.191 
0.  189 
0 . 189 
0.  189 
0.  189 
0.  188 
0.  187 
0.185 
0. 184 
0.182 
0.  177 
0.173 
0.  169 
0.  167 
0 . 162 
0.137 


. 500 
. 537 
. 642 
.856 
.067 
. 215 
. 281 
. 283 
. 246 
. 184 
. 118 
.983 
. 866 
771 
.693 
.632 
. 524 
.454 
.410 
. 379 
339 
. 318 
. 30  7 
. 299 
295 
. 294 
i.  293 
294 
.297 
.299 
i.  30  2 
i.  306 
i.  311 
i.  314 
I.  315 
.317 

1.317 

1.318 
i.  315 
I.  313 
1.312 
i.  308 
I 280 


0. 500 
0.531 
0.611 
0.788 

0. 969 
1.096 

1.  157 
1.161 
1.131 
1.080 
1.022 
0.902 
0.797 
0.710 
0.641 
0.584 
0.484 
0.420 
0.376 
0.347 
0 . 309 
0. 288 
0.276 
0. 268 
0. 265 
0. 262 
0. 261 
0. 260 
0.261 
0. 262 
0.265 
0. 268 
0. 270 
0. 271 
0. 272 
0. 271 
0. 271 
0. 269 
0. 267 
0. 264 
0. 261 
0. 258 
0. 23C 


0. 500 
0.525 
0 . 58  S 
0. 742 
0.900 
1.010 
1.063 
1.066 
1.040 
0.994 
0.942 
0.836 
0. 740 
0.661 
0.  597 
0.  545 
0.453 
0. 393 
0. 353 
0. 324 
0. 289 
0. 267 
0. 254 
0. 246 
0. 241 
0. 238 
0.  235 
0.  234 
0. 235 
0. 235 
0.  235 
0.237 
0. 239 
0. 240 
0. 240 
0.  240 
0. 239 
0. 237 
0 . 234 
0.  231 
0. 227 
0.  225 
0. 200 


0.500 
0.524 
0.578 
0.713 
0.852 
0.950 
0.996 
0.999 
0.975 
0.933 
0.886 
0.787 
0.699 
0.625 
0.565 
0.515 
0.426 
0.369 
0.329 
0. 30  3 
0.267 
0.248 
0.235 
0. 228 
0.  222 
0.  220 
0. 215 
0.214 
0. 214 
0.214 
0.  214 
0. 214 
0.214 
0.213 
0.212 
0. 211 
0 . 20  9 
0.  206 
0. 203 
0 . .'iOO 
0.  196 
0.  193 
0.  168 
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APPENDIX  I.  HIGHER  APPROXIMATIONS  TO  THE  TRANSPORT  COEFFICIENTS 


The  quantity  f _ appearing  in  Eq.  (6)  is  given  by  the  series 


f - =1+S  +T  + 

^ 1 2. 


-where 


bM  b22 


* * 


b„  (b;Sjb23  - b,3  ) 


( bu  • b"j  bu  b22  b33  + 'b,2b!3ba3  ~b«<  b*»  " ' h 33  ^ ) 


bu  - 4-  a 


(2,2)  * 


b = 7 n - 8 Cl 

12  ) 


bl3-f 

zz  * *■ 


K * n 321  (2,5)*  i35nU«*  ,.n(5,5)lt 

b23  IT  u -7-  0.  + TX1  ■ 30X1  , 

b33  = ^ 


The  quantity  f ^ appearing  in  Eq.  (7)  is  given  by  the  series 


f,  = l + 6,  + S.  + — , 

A 1 ^ 
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where 


'/  z 
a 

12 

II  U " « 

a a,,  - & 

ll  XX.  12. 


and 


//  . >/  « 'i  '/  \ - 

a (a  a -a  , a ) 

nv  ia  io 


( % -sv  **«  *x a » 4 <•  k *s.  - <;  c - ^ ^ - *;» «) 


a " = b , a"  = b , a"  = b , 

11  11  12  12  13  13 


nn  (1  })»  (a,3)  * _ (3*0  * 

" - n a -2sn  + 2on  , 

m — il  * u ' 


22  " 9 


" _ Sis  n(a,a5,t~  s<’1  nca-3)^  — d!9,4)*  - 30 n 


(3,5)#- 


„ I4S53  ^)*  ;aiff  n<aJ>*  + '545  ^ * ,os  q^**^*0* 


*33 


35& 


■n 


-rr 


The  third  approximation  to  the  quantity  f appearing  in  Eq.  (8) 

is  expressed  in  the  determinant  notation  of  Chapman  and  Cowling  as 

(2) 


f 


(3)  _ 


= a 


jj4 

00  J 


00- 


w 


(37) 


where  jtj  is  a fifth  order  determinant  whose  general  element  is  a-  , 
with  i and  j ranging  from  -2  to  2.  The  tern&n. is  the  minor  of 
M ^ obtained  by  deleting  the  row  and  column  containing  agQ  The 
elements  axj  are  functions  of  the  collision  integrals,  the  molecular  weights, 
and  the  'etative  con-,  er.tr ations  of  the  gases  in  the  diffusing  mixture  In  general, 


a.j  - r.^  ; and  if  i . J * 0 


a • has  the  same  form  as  a-  except  that  the 

J J 
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subscripts  1 and  2 , referring  to  the  molecular  species,  are  interchanged 


in  the  expression  for  a„  . The  a^  necessary  for  the  evaluation  of  r,q.  (37) 

are  as  follows,  where  x,  and  x_  are  the  mole  fractions,  and  M and 

12  12 

the  molecular  weights  of  species  1 and  2 , respectively: 

(!,')  * 


a =8 


M,Ma 


i n 


°o  12 


W*  ^ r 
a0i  8(M,.Ma)f 


5 - 3 o,3)* 

Tn„ 


a - ..  (»"■)* 

02  (M,  + Ma)^ 


35  _ 21  <■,»)*  , 

g 2 + 6'°')a 


o-i 


3 

iMj  aoi 


0-2 


(?;) 


a 


a. 


JS 

x. 


" ) 


(38) 


8M,  [ 5 . , a.  is  2 Ci,a)*  a ca.,3)* 


" (M..M/1  s 


" [(r„) 


-.a, 


I 
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a : Cl  + CL.  2 , 

a k.  12  - 


a 


*3  M . 


ia 


(W 


f ( a hN5«X’*'  1 <C 


X <l«o*  ‘2,2)*  „ „0k,i)* 

-IS  Ma,Q/’  +7M,M2Aa  -2K,Mani2 


(38) 
(cont. ) 


a 


2Ma  \a 


11  ' M,  + 


(2,2)*  „ (2.3>* 

7a  -8Xi, 


(O 


L-  rn)  J 
a 


8 ( M,Ma)  " 

a'-!" 


55  0,0*  16  «,2)*  0,«*  <-■»*  ' . 

TA»  ■ 2 “a  *4n.*  '2a'»  J ’ 


a - 

i- 


£ 2 
6M,S  M,1 


1 (M.4-W.) 

* .6 


f £95  cm)*  _ 547  a<'.*>*+  5£  /5a0>^ 

32.  12  IC.  iz  a 12,  I 


la 


.w™*  t8n';‘'3,* 


i a. 


12 


a - a.'  ■+  — cl 

22.  22  Xa  22-  ’ 


HM, 


r 


22 


(H,+  MJ 


(to  M(  t M M,X  ♦ 35  hX'  ^(«^35M;)n, 

12? 


0,2)* 

IX 


^M7i0«M>l33M><;f-l05<rff7t5HX'5>' 


12 


2 (2  2)*  3 _ 0.3)*  3 O,*-) * 

+•  ^M,na(4M1  +7Ha)Uja  -56M,M3vQ)a 


i a -,(3.3)* 
+ 12.  M M1AiJ 


J 


L"  s / 

32  if VV  11 


(rm> 


(OJ 
m a 


a. 


3(M,Mt)a 
:-2‘  ( m,+  K2)5 


(,4)*  cO0,£)* 


f 8505-0,1)*  2499  n(,>2)*  + — Xl'^*-  105  fL  ’ * a5^ 
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The  a.' . above  are  the  same  as  those  defined  in  connection  -with  the  thermal 

ij 

conductivity,  but  now  it  is  necessary  to  specify,  by  subscripts  on  the 

( J . n)* 

, which  molecular  interaction  is  under  consideration;  for 
example,  a^  refers  to  an  interaction  between  two  molecules  of  species  1 , 
whereas  a."  ^ ^ refers  to  species  2 . The  above  expression  for  the  a„  were 
calculated  by  the  method  of  Chapman  and  Cowling^^)  . 

The  expression  for  f^  may  be  considerably  simplified  for  the 
cases  of  tracer  diffusion  and  self- diffusion.  For  tracer  diffusion  where 
Xj — * 0 , we  obtain 


1 


(39) 


where 


01 


a * a 
Ma(tC  - 5) 


a a 
00  11 


30  M*  + 25M*-  f2n^e>*+lfeM,MaA* 


(a0i  ai£  ~ ^02.  k | | ) 
aooail^a(l  a22  ' a 12  ) 


The  quantities  A*  , B*  , and  C*  are  defined  in  Eq.  (11).  The  result  for 


x^  — * 0 is  obtained  from  the  above  by  changing  the  signs  of  all  subscripts  on 
the  a. ; , or  what  has  the  same  effect,  interchanging  M and  M-  . The 

1 w 


expression  for  self^diffusion  is  obtained  by  setting 
(In')# 

SI  ’ =fL.  , and  (r  ) = (r  ) = (r  ) . 

ii  22.  m 1Z  mil  m Lc 


si 


it,n)  * 


l 2 


(34) 


Reference  (4),  Chapter  9. 
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The  general  expression  for  the  thermal  diffusion  ratio  given  by 
Eq.  (9)  involves  the  a.,  defined  above.  For  a binary  mixture  of  heavy 

Ur 

isotopes,  the  dimensionless  quantity  k^,  , defined  in  Eq.  (10)  is  given 
to  a first  approximation  by 


o-i 


5 Cl  +2Cl  | 


v a a „ 

\ 1 1 ii 


n 


IS  (6C*  -5)  ( 3A*  + 5) 

i A*  (16  A*  - 12  & * + 5 5 ) 


(40) 


and  to  z.  second  approximation  by 


f Lx.-x.y.-x.yJ  , 


where 


Xl  = 


X2  “ 


3 Q-o-i  a a 5a~ 


0-2  aia. 


n 2 


_ n _ n 
Cl  Q.  - CL 
H 2*  1 2 


ao-i  ~ ^o-a,  aia 


a,i  Si 


i a 


X3  = 


^o-a.  it  ~ Vi  ^ i2 
/ . ' _/ 2 
aii  Si  ' aia 


Y,  = 


c,  aaa  - cA  aa 


n - n - //  3 

% a32  ' ttii 


YZ  = 


^3  ^ ~C3ai2 

« u n z 

an  a22’  aia 


1 


? 


3 


(41) 
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and 


45  0,i)  * 4.5  (is)*-  _ 

C = 7 A - — A t 18,0. 
1 ? 


C-,  = 


1445  ^ l? 2.7  225  ^(l’3)  " 75 


C3  ' 


512 


P\ 

. A L 


Cl,^-)  X- 


2 " 64  32  - i - 2 

= 7805  no,'>*  5*Mxitl,a)>  + S3iq  j0°,3)*  5affnci,‘°*  i i£f  a(i,5)* 


64 


16 


2C  (2,2)*  (3,3)*-  (3  +)*- 

_ 35  ^ ■ _ |+ri  + ion  +30, 


1+0 


In  the  above  expressions  for  k^,  the  are  to  be  calculated  as  if 

Mj  = 1 inasmuch  as  the  deviations  from  this  equality  have  been  explicitly 

allowed  for  in  deriving  Eq.  (10)  . It  is  interesting  to  note  that  for  rigid 
spheres  [k^,]  ^ Is  8.  8 percent  larger  than  [k^,]  ^ . 

The  above  expressions  for  the  higher  approximations  to  the 
transport  coefficients  were  calculated  by  the  method  of  Chapman  and  Cowling^®), 
and  their  form  depends  on  a particular  approximation  scheme  for  obtaining 
solutions  to  an  infinite  set  of  simultaneous  equations.  The  scheme  used  by 
Chapman  and  Cowling  is  not  unique,  and  indeed  Kihara^^)  has  developed  an 
alternative  scheme  and  from  it  obtained  the  following  second  approximations 
to  viscosity,  thermal  conductivity  r and  self -diffusion,  and  the  first 


approximation  to  the  thermal  diffusion  ratio  for  isotopic  mixtures; 
.(a) 


. (2,3)*  - ^ 


f ~,  1,  -i 

^ fq  l a12'2'*  a.  J 1 


(42) 
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;s) 

C = 

Ji 


4 XL 
si 


(5,5)  * 


(43) 


«>  ( , q*  _ 5 ) 

p v u>  ^ ->  y 

V * | + 

D It  A*  + *r0 


(44) 


15  (fc  C.*  - 5 ) 

Mo  A* 


(45) 


Kiiuira'*  expressions  are  considerably  simpler  than  the  corresponding 
Chapman  and  Cowling  expressions,  and  are  known  to  be  more  accurate  for 
the  special  case  of  a Lorentzian  gas  (a  binary  mixture  in  which  one  component 
has  a much  smaller  molecular  weight  than  the  other,  and  collisions  between 
molecules  of  the  light  component  need  not  be  considered).  It  is,  therefore, 
of  interest  to  compare  the  expressions  of  Kihara  and  of  Chapman  and  Cowling 

I 

for  the  more  realistic  potential  given  by  Eq.  (1),  and  to  estimate  their 
relative  accuracy  by  comparison  with  the  still  higher  approximations 

j 

evaluated  in  the  present  paper. 

Reference  (4),  Chapter  7. 

T.  Kihara,  Imperfect  Gases  (originally  published  in  Japanese  in  1949, 
and  translated  into  English  by  the  U.  S.  Air  Force),  Secs.  23,  26,  27, 

31,  32;  "Virial  Coefficients  and  Models  of  Molecules  in  Gases", 
University  of  Wisconsin  Report  OOR- 7,  June  5,  1953,  to  be  submitted 
to  Rev.  Mod.  Pbys.  for  publication. 
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The  two  methods  yield  answers  for  the  potential  o£  Eq.  (1)  which 
are  appreciably  different  only  for  the  case  of  the  thermal  diffusion  ratio; 
some  selected  values  for  isotopes  are  given  in  Table  EX  below.  The  overall 
accuracy  of  Kihara's  expression  seems  to  be  somewhat  better.  In  view  of 
this  and  of  their  greater  simplicity  (fewer  different  collision  integrals  are 
required),  Kihara's  expressions  for  the  second  approximation  to  viscosity 
and  thermal  conductivity  and  the  first  approximation  to  the  thermal  diffusion 
r#tio  would  seem  to  be  definitely  preferable  to  Chapman  and  Cowling's 
corresponding  expressions.  It  is  for  this  reason  that  we  have  avoided  expressing 
our  results  for  isotopic  thermal  diffusion  in  terms  of  the  usual  quantity  R^  , 
defined  as  the  ratio  of  to  the  first  approximation  to  for  rigid  spheres. 

Since  the  first  approximations  to  k,^  for  rigid  spheres  according  to  Chapman  and 
Cowling  and  to  Kihara  are  different,  R^,  involves  an  uncertain  numerical 
factor,  and  the  quantity  k^.  defined  in  Eq.  (10)  seems  preferable  from  a 


theoretical  standpoint. 


19  June 


Table  IX  Comparison  of  the  approximations  of  Chapman  and  Cowling 


PHH 


141  2 

CC  K 

CC 

0. 3062 

0. 3104 

0.3133 

• 0.0828 

-0.0819 

-0.0829 

•0.0465 

-0.0469 

-0.0456 

0. 1860 

0. 1903 

0. 1841 

0 . 4349 

0. 4455 

0.4536 

0.4197 

0.4291 

0.4387 

0.3852 

0. 3942 

0.4024 

OL  = 15 

0. 3062 

0. 3104 

0.  3133 

0.0217 

0.0216 

0.0229 

0.0505 

0.0511 

0.0498 

0.2668  ‘ 

0. 2739 

0 . 2665 

0. 5189 

0. 5348 

0.  5446 

0. 5204 

0. 5357 

0 . 5479 

0.5253 

0. 5422 

0. 5534 
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APPENDDC  II.  COMPARISON  WITH  AMD  UR'S  SEMI-CLASSICAL 
CALCULATION  OF  S(  ^ ^(K) 

( 3 T) 

Amdur  has  devised  an  approximate  method'  for  the  calculation 
of  gaseous  transport  properties  which  is  capable  of  handling  very  complicated 
inter  molecular  potential  functions  with  relatively  little  computation,  and  has 
applied  his  method  to  the  calculation  of  the  transport  properties  of  the  rare 
gases  at  low  temperatures.  The  procedure,  however,  involves  assumptions 
which  require  further  justification^3®^  , so  that  it  is  of  interest  to  provide  a 
test  of  the  assumptions  by  comparing  the  results  of  the  present  calculations 
with  those  obtained  by  applying  Amdur ‘s  method  to  the  same  intermolecular 
potential. 

In  the  approximate  procedure  the  transport  cross  sections  are  not 
calculated  directly,  but  instead  the  total  classical  collision  cross  section  is 

first  obtained,  and  the  transport  cross  sections  are  then  assumed  to  be  in  the 

» 

same  ratio  to  the  total  collision  cross  section  as  in  the  case  for  rigid  spheres. 

The  total  collision  cross  section  is  infinite  if  calculated  by  completely  classical 

means,  so  a semi-classical  approximation  is  used,  in  which  only  scattering 

through  angles  greater  than  a certain  critical  angle,  X , is  considered.  The 

c 

*37)  I.  Amdur,  J.  Chem.  Phys.  _15,  482  (1947);  Jh,  190  (1948);  JJ,  100 
(1949). 

^3®)  H.  S.  W.  Massey  and  E.  H.  S.  Burhop,  Electronic  and  Ionic  Impact 

Phenomena  (Oxford  University  Press,  London,  England,  195>2).  pp.  381- 
382,  385. 


total  collision  cross  section,  S(K)  , is  then  given  by 
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S(K)  = 7i  b*  , (46) 

■where  b_  is  that  value  of  the  collision  parameter  which,  for  a given  value  of 
c * 

the  energy  K , will  produce  a deflection  equal  to  X • The  velocity 

c 

dependent  cross  sections  for  viscosity  and  diffusion  are  then  taken  to  be  two- 
thirds  and  one-half,  respectively,  of  S(K)  , as  for  rigid  spheres.  The  collision 
integrals  are  obtained  from  the  cross  sections  by  essentially  the  same  procedure 
as  used  in  the  present  paper. 

Amdur's  method  was  applied  to  the  potential  given  in  Eq.  (1)  to 

calculate  approximate  values  of  the  velocity  dependent  transport  cross  sections 

for  comparison  with  the  accurate  values  obtained  in  this  paper  by  direct  means. 

(39) 

The  critical  scattering  angle  was  taken  to  be'  ' 

v . _A__. 

c 2yK*  5 (47) 

where  A is  the  dimensionless  group  h/  [ rm(2/t  £ ] » h is  Planck's 

constant,  and  the  other  quantities  are  as  previously  defined.  Following  Amdur, 

X "was  represented  by  Eq.  (16)  with  the  second  term  neglected;  use  of  this 
c 

equation  restricts  the  range  of  the  results  to  K > 1 . The  results  are  expressed 
in  reduced  form,  S*(K)  = S(K)  /(  fr  r*  ) . According  to  the  approximation, 

S*(K)  = ^ \yL)  for  all  values  of  JL  : however  for  K > 1 , is  is  actually  found 


(39) 


H S.  W.  Massey  and  C . B.  O.  Mohr,  Proc.  RoySoc.  (London)  A141, 
434  (1933). 
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that  S^(K)  is  always  larger  than  S^'iK)  , by  as  much  as  49  percent. 

In  Fig.  5,  S*(K)  is  shown  for  f\  =0.05  and  2 , where  the  potential 

parameter  oL  = 12  The  small  value  of  A corresponds  to  a heavy  gas  -with 

strong  inter  mo  leeula  r attraction;  such  as  xenon,  the  large  value  to  a light  gas 

with  weak  attraction,  such  as  hydrogen  Also  shown  for  comparison  are 

S^(K)  and  S^(K)  for  oc  - 12.  The  results  for  other  values  of  oc  are  quite 

similar.  It  is  seen  that  although  S*(K)  does  not  change  rapidly  enough  with 

K , it  maybe  nearly  correct  over  a small  range  of  K . This  fact  would 

account  for  Amdur's  numerical  results  on  the  rare  gases  (except  helium), 

where  he  found  fair  agreement  between  calculated  and  experimental  values  of 

the  viscosity,  but  the  calculated  temperature  variation  was  too  small.  The 

calculated  values  of  the  self-diffusion  coefficients  of  neon  and  argon  were, 

however,  rather  too  small^®^^,  as  might  be  expected  since  ^(K)  is 

really  smaller  than  S^(K)  . The  agreement  for  helium  was  somewhat  better, 

but  this  was  perhaps  due  in  part  to  the  use  of  a potential  which  is  no  longer 

(42) 

considered  very  accuiate'  ' 

The  approximation  scheme  is  thus  seen  to  give  results  which  are 
valid  over  a limited  range,  and  hence  may  still  .prove  useful  because  of  its 
ability  to  handle  very  complicated  potential  functions, 


(40) 

(41) 

(42) 


E.  B.  Winn,  Phys.  Rev.  8(3,  1024  (1950). 

F.  Hutchinson,  J.  Chem.  Phys,  1 7,  1081  (1949). 
Reference  (36),  p.  378. 


Jig,.  5.  Comparison  of  exact  and  approximate  calculation*  of  the  reduced,  crosa  sections, 
S(K),  for  the  Buckingham  potential.  Solid  lines:  calculated  ly  numerical  integration; 

broken  lines:  calculated  by  Amdur's  approximation  scheme. 


